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Abstract. We prove the existence of minimizers for Hartree-Fock-Bogoliubov 
(HFB) energy functionals with attractive two-body interactions given by New- 
tonian gravity. This class of HFB functionals serves as model problem for self- 
gravitating relativistic Fermi systems, which are found in neutron stars and 
white dwarfs. Furthermore, we derive some fundamental properties of HFB 
minimizers such as a decay estimate for the minimizing density. 

A decisive feature of the HFB model in gravitational physics is its failure 
of weak lower semicontinuity. This fact essentially complicates the analysis 
compared to the well-studied Hartree-Fock theories in atomic physics. 



Contents 

1. Introduction 

2. Basic Properties of HFB Energy 

3. Main Results 

4. Preliminaries 

5. Proof of Theorem [2] 

6. Proof of Theorem [3] 

7. Proof of Theorem [1] 

8. Proof of Theorem H 

Appendix A. Localization of Kinetic Energy 
Appendix B. Proof of Proposition [2TT] 
Appendix C. Proof of Proposition |42] 
References 



1 

6 
9 

14 
2C 
27 
28 
34 
36 
38 
39 
42 



1. Introduction 

The Hartree-Fock-Bogoliubov (HFB) theory is a widely used tool [37l[T2] for un- 
derstanding many-body quantum systems where attractive two-body interactions 
are dominant. In particular, HFB energy functionals incorporate the physical phe- 
nomenon of "Cooper pairing" which is likely to occur whenever attractive forces 
among quantum particles become significant; e.g., in nuclear and gravitational 
physics as well as in superconducting materials. Despite the broad range of phys- 
ical applications for HFB theory, not much has been known rigorously concerning 
the existence of minimizers, let alone a proof of their fundamental properties such 
as Cooper pair formation. 
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As a starting point for rigorous analysis, this paper is devoted to the existence 
of niinimizers for HFB relativistic energy functionals with an interaction which 
behaves at infinity hke the (attractive) Newtonian gravity. This specific class of 
functionals can be viewed as a model problem for self-gravitating relativistic Fermi 
systems which are found, for example, in neutron stars and white dwarfs. Moreover, 
speaking from a mathematical point of view, the corresponding variational problem 
exhibits the delicate property of criticality, as we will detail below. 

The most challenging main feature of the HFB variational problem in gravita- 
tional physics is its lack of weak lower semicontinuity (wise) due to the attractive 
interaction among particles. As a consequence of the absence of wise, the existence 
proof for minimizers is much more involved than for the well-studied Hartree-Fock 
(HF) models arising in atomic physics, where wise plays an essential role; see 
[30l l36] . Another difiiculty in the analysis of HFB models stems from its trans- 
lational invariance and the fact that the main variables are two operators related 
via a complicated constraint inequality. Finally, a further complication (although 
conceptually less important) is the treatment of the pseudo-differential operator 
describing the kinetic energy of relativistic fermions. 

Following the general rigorous discussion of HFB energy functionals in [2], we 
specifically consider the energy functional given by 

(1.1) f(7,a) :=Tr(T7) + | //' W{x - y)p^ix)p^{y)dx dy 

-J// W{x-y)\j{x,y)\^dxdy+'^ W{x - y)\a{x,y)\^dx dy, 

which models a system of relativistic fermions subject to the interaction W. We will 
assume that W{x) behaves like — l/|a;| at infinity, i.e. that far away the attractive 
Newtonian force is predominaniQ. Indeed in most of the paper we will even take 
for simplicity 

In (jl.ip . Tr(-) denotes the trace, and the pseudo-differential operator 



T — v— A + m^ — m 

describes the kinetic energy of a fermion with mass m > 0. The coupling constant 
K > parametrizes the strength of the interaction among the particles. Further- 
more, the variables 7 and a are two operators acting on L^(R'^; C^)0 The operator 
7 is referred to as the one-body density matrix. It is self-adjoint, nonnegative and 
has a finite trace fixed to be 

(1.2) Tr(7) = A, 

for some given A > 0, which is a real parameter and it corresponds to the expected 
value of the number of particles in the star. The density p^ is the unique L^(K'^) 
non-negative function that satisfies TT{'-fV) = p-y{x)V{x) dx for any bounded 
function V. As 7 is trace-class, it has a kernel j{x,y) appearing in the second line 
of (jl.ip . which is a 2 X 2 hermitian matrix for almost every x,y eM.^. Finally, the 



^See Remark [6] for precise assumptions on W. 

^More generally, we shall consider L^(R'^;C) below, where q ^ 1 denotes the internal spin 
degree of freedom. The case q = 2 (corresponding to spin 1/2) is physically the most relevant one. 
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operator a is called the pairing density matrix. It is only assumed to be Hilbert- 
Schmidt, that is we have TT{a*a) < oo. Its kernel is a 2 x 2 matrix which is 
supposed to be antisymmetric in the following sense: a(x,y)'^ = —a{y,x), where 
is the usual transposition of matrices. 

To formulate the HFB variational problem for the energy £(7,0;), we have to 
supplement the side condition (|1.2p by the following operator inequality relating 7 
and a: 

<i-^" (0 o)«(:- 1-7) ;) »i^r;c')eL»(M-c"). 

This inequality guarantees that the pair (7,0) is associated to a unique quasi- free 
state in Fock space [5]. The corresponding HFB minimization problem then reads 

J(A) := inf {£:(7,a) : 7* = 7, = ~a such that dT!]) and dO]) }. 

For the physical interpretation and some background from many-body quantum 
mechanics, we refer the reader to [3J[3]. 

Our first main result is formulated in Theorem [T] below, which shows that 

/(A) is attained for A < A"^^(k) and k < 4/7r. 

More precisely, we prove that all minimizing sequences (with such A and k) are 
relatively compact, up to translations, in an appropriate topology for the pair of 
operators (7,0;). In physical terms, the finite number A^^^(k) can be interpreted 
as the Chandrasekhar limiting mass for the HFB model considered here. That is, 
a self-gravitating relativistic Fermi system with supercritical particle number, 

A> A«F^(.«), 

can undergo "variational collapse" in the sense that £(7, a) fails to be bounded 
below and hence minimizers cannot exist in this case. Regarding the relation to 
semi-classical models of stars, we remark that a useful estimate for A^^^(k) is 
provided by the asymptotic estimate (when the interaction is purely gravitational): 

Here Tc — 0(1) is a numerical constant that can be calculated from the classical 
Thomas- Fermi-type theory for neutron stars and white dwarfs, which was initiated 
by S. Chandrasekhar in his seminal work [5]; see also [?T1 132) . 

Apart from the existence of minimizers, we also establish some of their fundamen- 
tal properties. In particular, we derive a certain decay estimate for the minimizing 
density p-yix), which turns out to be essential when proving the existence of mini- 
mizers, as we will detail below. Here, we just point out that the fact that a decay 
estimate for pj cannot be simply obtained by using the self-consistent equation for 
HFB minimizers, since the minimizing one-particle matrix 7 is not known to have 
a finite rank. In fact, by using results from [1] [T3], we shall see below that any 
minimizer (7, a) in the spin-1/2 case exhibits an infinite rank, provided we assume 
a 7^ holds. To deal with this difficulty arising from infinite-rank properties of 
HFB minimizers, we devise some variational arguments combined with a second 
order expansion of the HFB energy to prove an estimate for the fall-off behavior of 
the minimizing density p-y{x). 
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Concerning further fundamental properties of minimizers for £(7,0), we remark 
that it is an important open question to show that minimizers actuahy exhibit 
non-vanishing pairing a 7^ 0, at least for a coupling constant k which is not too 
small (relatively to A~^^'^). On heuristic grounds, one expect such a phenome- 
non of "Cooper pair formation" to be energetically favorable due to the attractive 
interaction among particles. However, it seems to be a formidable task to find 
mathematical proof for this claim. 

Let us briefly comment on our existence proof. The relative compactness of 
minimizing sequences (and, in particular, the existence of minimizers) is the con- 
sequence of the validity of the binding inequality: 

(1.4) /(A) < /(A') + /(A - A'), for all A' e (0, A), 

which is also referred to as strict sub-additvity condition. In fact, binding inequali- 
ties like (|1.4|) appear naturally in the analysis of compactness properties of minimiz- 
ing sequences; for instance, when using the "concentration-compactness principle" 
as done in [34l[35]. Moreover, binding inequalities also turn out to be useful for 
linear models, in which the bottom of the essential spectrum has the form of the 
minimum on the right side in (|1.4[) , as expressed by the celebrated HVZ Theorem; 
see [23l|43l|44l|13l|40l[39]. The interpretation of is that it is not favorable for 
a minimizing sequence to have a mass A — A' escaping to infinity, while the mass A' 
stays in a bounded set (up to translations). 

However, there is a notable difficulty if one tries to prove (ll.4p for the HFB model 
studied in this paper. To see this, we recall that in the usual TV-body Schrodinger 
setting, the binding inequalities are quantized, i. e., they are of the form 

/sch(iV) < /sch(iV - N') + /sch(iV'), for all iV' = 1 . . . - 1. 

This quantization just refiects the fact that we can restrict ourselves to states with 
a definite particle number and use adapted geometrical methods [13l [40l [39] . In 
particular, an important feature of quantized binding inequalities is that they can be 
proved by induction; see, e.g., [TS]. In HFB theory, by contrast, it is not expected 
that ground states will have a definite particle number. Therefore, the binding 
inequality ()1.4p cannot be simply established by induction in this model. 

To overcome this difficulty and eventually prove the binding inequality ()1.4p for 
the HFB model considered here, we use a different method in the spirit of a tech- 
nique introduced by P.-L. Lions in |36| . More precisely, we assume the existence of 
a non-relatively compact minimizing sequence with particle number A < A^^^(k), 
and we describe its behavior in detail. In particular, we show its splitting into 
several pieces receding from each other, where at least two parts are relatively com- 
pact and carry strictly positive particle numbers A^ > and A^ > 0, respectively. 
Furthermore, we conclude that the infima /(A*) with i = 1, 2 must be attained, and 
we find that the ground state energy decomposes as 

/(A) = /(Ai) + /(A2) + /(A - Ai - A2). 

On the other hand, by deriving a suitable decay estimate for any minimizer of /(A'), 
we can deduce the strict inequality 

/(A) < /(A^) + /(A^) + /(A - A^ - A^), 
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in contradiction to the previous equality. Thus the main assumption (failure of 
relative compactness of all minimizing sequences) cannot hold. This rationale es- 
tablishes a posteriori the binding inequality (jl.4p for the HFB model. 

Although the aforementioned behavior of minimizing sequences appears phys- 
ically quite clear, implementing the above argument mathematically is far from 
being an easy task. One difficulty stems from the fact that the main variable 
is a pair of two operators satisfying the complicated constraint (|1.3|) . whereas in 
usual nonlinear PDE problems, the variable is often a single scalar or vector-valued 
function. In order to deal with this issue, we use some ideas of previous works 
[l8l [T9l [T4] . Also, we point out that the pairing term a cannot be handled by 
obvious means: In contrast to the trace-class operator 7, the operator a is a priori 
only Hilbert-Schmidt. Moreover, loosely speaking, the direct and exchange term 
are both "suhcritical" in the sense that they may be controlled by the L^^^^-norm 
of ^ e H^/'^{9?). By contrast, the pairing energy depending on a is "critical" 
because it can only be controlled by the kinetic energy of 7 itself. A particular 
illustration of this difficulty is that even ruling out the vanishing of a minimizing 
sequence is quite delicate, and that the proof of the decay of minimizers is based 
on two-body techniques for the wavefunction a{x, y). 

Organization of the Paper. The outline of this paper is as follows. Section[2]sets 
the stage for the variational calculus related to the HFB energy (7, a) i~?> £(7, a). 
Our main results are then formulated in Section[3] After collecting some preliminary 
facts in Section |4l the proofs of the main theorems are presented in Sections [SHE] 
below. The appendix contains various technical results and proofs. 

Notation. We employ standard notation for and Sobolev spaces. For an integer 
g ^ 1, we define the inner product (/, 5) = /^a Ylt=i fii^)9ii^) dx for functions / 
and g in ^^(R^; C). If M is a g x g matrix, we use the notation jMp = Tr(M*M). 
In the whole paper we employ the physicists' notation for the operator 

/ I— > {Lp,f)ip. Furthermore, we shall use 6p, with 1 ^ p < 00, to denote the 
Schatten class of operators A acting on i^(R'^; C) and having a finite p-trace, and 
we endow it the norm ||A||g Tr(|^|P)i/P < 00; see gT]. The space ©00 denotes 
the space of compact operators on L^(R^;C'^), equipped with the operator norm 
II • ||. Below we shall also introduce a "Sobolev-type" Schatten space A" c ©i x 62 
in the study of the HFB energy functional 8. 

For later convenience, we introduce two smooth cutoff functions XR Cr that 
are defined as follows. Let ^ x ^ 1 be a fixed smooth function on R'^ such that 
X = 1 for |a;| < 1 and x = fo'' N| ^ 2. For any i? > 0, we then define the functions 



Throughout this paper, we use the cutoff functions XR ^nd C_r freely. Furthermore, 
we denote by Ty : L^(R^; C^) — > L^(R^; C) the unitary operator that is defined by 



for given y € R^. We also use the notation U{"f,a)U* :— {UjU* ,UaU*) for any 
unitary operator U acting on L^(R'^; C), e. g., for the translation operator U ~ Ty. 

For the physically inclined reader, we remark that we work in units such that 
Planck's constant H and the speed of light c satisfy h — c — 1. 



(1.5) 



Xr{x) = x{x/R) and Cr{'x) = Vl - Xr{x)'^- 



(1.6) 



{ryf) = f{-y), 
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2. Basic Properties of HFB Energy 

To prepare the statement of our main results, we first provide the adequate 
setting for the Hartree-Fock-Bogoliubov (HFB) variational problem, and we collect 
some basic properties needed for the rest of this paper. 

2.1. HFB States With Finite Kinetic Energy. To set up the variational calcu- 
lus, we define a class of HFB states having finite pseudo-relativistic kinetic energy. 
Therefore we introduce the following (real) Banach space of density matrices 

(2.1) X = {(7,a) e ei X 62 : 7* = 7, = -a, ||(7,a)lk < 00}, 
equipped with the norm 

(2.2) ||(7,«)IU = 11(1 - A)i/S(l - A)i/ig^ + 11(1 - A)i/4a||g^. 

We remind the reader that 61 and 62 denote the space of trace-class and Hilbert- 
Schmidt operators on L^(IR.^; C), respectively. The universally fixed integer g ^ 1 
takes into account the internal spin degrees of freedom of the model. In what 
follows, we shall omit the dependence on q whenever it is of no importance. 
Furthermore, we define the following subsets of density matrices in X: 

(") (2 »)<(:: ,_%).(; \)\ 

and, for A > given, we set 

(2,4, (° °) . ?) , IVh) ^ a} , 

Here 7 is defined by complex conjugation of the kernel j{x,y). Note that both 
JC and IC\ are closed and convex sets in X, and that we have ^ 7 ^ 1 for all 
(7, a) e /C. In the following, we shall use the notation: 

(2.5) Tr(r7) := Tr ((-A + m2)i/S(-A + m^/^^ - mTr(7), 

for T = A + — m. Also, we sometimes allow the case of vanishing mass 
m — 0. But we remind the reader that strict positivity of to > is an essential 
assumption in all the main theorems stated below. 

As a first simple fact, we obtain the following Sobolev estimate for the square 
root of the density function p^. 
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Lemma 2.1. For any (7, a) G JC, the density function p-y satisfies 

Moreover, the map (7, a) f-^ \/P^ continuous from K- to LP{M.^) for all 1 ^ p < 3. 
Proof of Lemma [K7[ Writing 7 — J2n ^nl'fin) {'Pnl with ^ A„ ^ 1, we find 

(2.6) Tr{V^l) = ^ Xni^n, V^^n). 

n 

Furthermore, we recall that, for any (p^ip G H^^'^{M.^), 

(2.7) {VW+W,^VW+W)^{^,^^) + {^,V^^^)- 

This follows from [551 Theorem 7.13] and {\(p\, ^/-A\(p\) {(p, ^/-Aip). These two 
estimates are easily seen to imply the inequality in Lemma 12.11 The rest simply 
follows from Sobolev's embedding. □ 

2.2. Boundedness From Below of HFB Energy. We define the (purely gravi- 
tational) HFB energy as 

(2.8) f (7,a) := Tr (r7) - ^Dip,,p,) + ^Ex(7) - J // ^-^^dxdy, 

2 ' ' 2 2 77js.3xr3 \x-y\ 

where we now use the shorthand notations 

(2.9) Dif,g):^[[ Ij^d. dy, Ex(7) // ^^^d. dy, 

which we refer to as the direct term and exchange term, respectively. The last term 
on the right side of (|2.8p is called the pairing term. 

For (7, a) G /C, we deduce from Lemma [2.11 and the Hardy-Littlewood-Sobolev 
inequality that 

(2.10) D{p.„p^)i^C\\^p^\\l^y,<^. 

Thus the direct term D{pj,py) is well-defined. Next, we recall the pointwise esti- 
mate 

(2.11) |7(x,y)|2 <;p^(x)p^(y), fora.e. (x,y)e]R3xR3^ 

which is obtained by writing the spectral decomposition of 7 and using the Cauchy- 
Schwarz inequality for sequences. This yields 

(2.12) [ [ ^"f^ '^'^j dxdy ^ D{p^,p-.) < 00, 
JJr3xR3 \x-y\ 

showing that the exchange term is also a well-defined quantity when (7, a) G /C. 
To deal with the pairing term in £(7, a), we recall the Hardy-Kato inequality 

(2^13) R*f^- 

see [24l[22]. Applying ()2.13p in the variable x with y fixed, we obtain 
(2.14) / / \^^hy}}^dxdy < ;Tr(V^aa*) < 00, 

JJr^xr^ \x-y\ 2 ^ 

for any (7, a) G K.. Next, we introduce the 2 x 2-matrix 

■ 1-7 
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which defines an operator acting on L^(R'^;C'') L^(R'^;C''). Following the con- 
vention in [2 , we refer to any such F satisfying the constraint (|1.3p as an admis- 
sible 1-particle density matrix (1-pdm). In particular, we immediately obtain from 
^ F ^ 1 that F^ ^ F, and this leads to the operator inequality 

(2.15) 7^ +aa* 7, 

as already observed in [5]. This fact combined with (|2.14p yields 

k(x,2/)|2 



(2.16) ^^d.d,^5Tr(y3A7). 
Hence the pairing term is controlled by the kinetic energy. 

Remark 1. Because of (j2.16p . we will have to assume that ^ k < 4/7r m the fol- 
lowing; see also Remark\^ below why this upper bound on k is optimal. With regard 
to its physical application to stars, such a smallness condition for k is harmless, 
since we have k = Grn^ ^ 10~^^, when G is Newton's gravitational constant and 
m the neutron mass. Moreover, it is also natural to consider the HFB model in the 
limiting regime when k ^ 1 with kX"^/^ being fixed (i. e., the semi- classical regime 
for relativistic fermions). 

Remark 2. The arguments presented above can be easily adapted to prove that 
f : /C — 7- M is continuous with respect to the norm \\ ■ 

All this proves that the HFB energy functional (7,0) i— >■ £(7, a) with mass 
parameter m and coupling constant k is a well-defined functional on /C. Hence it 
makes sense to consider the following minimization problem 



(2.17) 



/(A) = inf {£(7, «): (7, a) G/Ca} 



with A ^ given. Note we have /(O) = with unique minimizer (7, a) = (0, 0). 

Before stating our main results on the variational problem (|2.17l) . we need to 
define the largest possible number A of particles of the system for which the energy 
is bounded below. We have the following statement. 

Proposition 2.1. (Chandrasekhar Limit). Let to ^ and ^ k < 4/tt be 
given. Then there exists a unique number A^^^(k) > 0, which is independent of m, 
such that the following holds. 

{i) For A s$ A"™(k), we have /(A) > -00. 
(ii) For A > A^^^(k), we have /(A) — —00. 

Furthermore, the function A^^^ (k) is nonincreasing and continuous with respect 
to K. It satisfies the asymptotic estimate 

(2.18) AH™(«:)^ as ^ 0, 

for some universal constant — 2.677. 

The proof of Proposition 12.11 is given in Appendix [B] below. The constant Tc 
was already defined in |32) . where it appears from a Thomas- Fermi- type energy 
obtained in the semi-classical limit k with nX^/^ — 0{1). 

Remark 3. It will be a consequence of our proof that /(A) = —00 for all X > 
and all m ^ when n > 4/7r; see Remark \15\ below. Saying differently, we have 
A"™(k) = when n > A/n. 
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Remark 4. One simplified model consists in restricting the variational set to HFB 
states without pairing, i. c, states of the form (7, 0) with ^ 7 ^ 1 and Tr(7) — X. 
A further simplification consists in restricting the energy to projectors only, i. e., 
7^ = 7, when X = N is an integer (this is the usual Hartree-Fock model). The 
so-obtained functionals satisfy similar properties as £(7,0:) on the corresponding 
variational sets, except that the assumption k < A/n can be removed. 

Remark 5. It can be easily seen from our proof provided in Appendix\B^ that for 
every ^ t < Tc 



where Ichi^, 1"^) is the Chandrasekhar ground state energy as defined in [3HI32j. 

3. Main Results 

In this section, we state our main results for the HFB variational problem defined 
in (|2.17p . First, we establish the existence of minimizers in Theorem [1] Essential 
properties of minimizers will then be formulated in Theorem [5] for the general case, 
and in Theorem [3] for the physically relevant q — 2 case. Furthermore, we study 
the simplified Hartree-Fock model and derive an existence result of minimizers; see 
Theorem m below. Finally, we briefly discuss the time-dependent HFB theory. 

3.1. Existence of Minimizers in HFB Theory. We first formulate the existence 
result of minimizers with particle numbers below the critical threshold. Since the 
functional £ is translation invariant, we have to take into account the unitary action 
Ty on L^(R'^) given by the group of translations in K'^; see Equation (|1.6p above for 
the definition of t. The precise existence result now reads as follows. 

Theorem 1. (Existence of Minimizers). Fix the integer q ^ 1 describing the 
internal spin degrees of freedom. Furthermore, suppose that m > and < k < 
4/7r. Then, for all < X < A^^^(k), the following properties hold. 

(i) Every minimizing sequence {(7n,an)}neN for /(A) is relatively compact in 
X up to translations. That is, there is a sequence {?/„} 

ngN C K'^ such that, 

after passing to a suitable subsequence, we have 

'''yS'^'"--' '^ri)Ty,^ — ?> (7, a) strongly in X as n ^ 00, 

where (7, a) G ICx is a minimizer for /(A). In particular, there exists a 
minimizer (7, a) for /(A). 
{ii) The following binding inequality holds for all < A' < A.' 

/(A) < /(A- A') + /(A'), 

Let us make some remarks. 

Remark 6. As can be checked from our proof, Theorem[l\ also holds true when the 
purely gravitational interaction —l/\x—y\ is replaced by a radial potential W{\x—y\) 
satisfying the following assumptions 

\W(\x\)\i^-^ forallxeR^. and VF(|a;|) sC --^ forlxl^Ro, 

\x\ ' 

for some constants e > and Rq > 0. Of course, the largest particle number 
A^^^(k) has to be defined accordingly. 



10 



E. LENZMANN AND M. LEWIN 



Remark 7. In fact, Properties (i) and (ii) are found to be equivalent, as it is 
well-known in the general framework of concentration- compactness method [34] for 
variational problems with translation invariance. 

Remark 8. Note that we could not prove the existence of minimizers for A = 
A^^^(k), for which the energy is nevertheless bounded from below. The reason 
is that there seems to be a lack of coercivity in this case (see Lemma \4-^ below). 
Furthermore, if m > Q holds, we do no expect minimizers to exist with critical mass 
A = A^^^(k). 5*66 also !32] for a related nonexistence result for a scalar Hartree 
equation describing boson stars. 

Remark 9. It is essential to assume that m > holds, which is also the physically 
relevant scenario. In the case of vanishing mass m — Q, a further loss of compact- 
ness (besides the one induced by translation invariance ) arises from scaling. 

3.2. Properties of HFB Minimizers. 

3.2.1. Nonlinear Equation and Decay Estimate. Let us now describe some proper- 
ties of HFB minimizers (7, a) obtained in Tlieorem[T] As before, we denote 



and we introduce ttie following HFB mean-field operator 
(3.2) Ft := 

acting on L'^{M?; C«) ® L'^{M? ;£''). Here 



\x-v\ 




(3.3) H, :^T- Kip, * ^)(x) + 

I- 1 \x~y\ 

is the usual mean-field operator of Hartree- Fock theory [30l |36l [2]. Moreover, it 
turns to be convenient to define 



(3.4) N = 



1 

-1 



Note that F commutes with N if and only if a = 0, that is if and only if the 
corresponding quasi-free state in the Fock space also commutes with the number 
operator J\f; see, e.g., [1] for the Fock space formalism in HFB theory. 

We can now state some fundamental properties of the minimizers for the HFB 
model considered in this paper. 

Theorem 2. (Properties of Minimizers). Let q ^ 1 be given and suppose 
TO > and < K < 4/7r. Assume that (7, a) G IC\ is a minimizer for /(A) for some 
A > 0. Then there exists a negative real number /i < such that T = r(7, o) solves 
the following nonlinear equation 

(3.5) T = X[~oofi){Fr- tiN) + D, 

where D is a finite rank operator of the same matrix form as T and satisfies 
ran(Z?) C ker(i^r — t^N). If moreover a 7^ holds, we have the improved bound 
fi < l3 — m, where the number (3 < m is given by Proposition^. S\ below. 
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Finally, the following decay estimate for the density function holds, for all R > 
sufficiently large, 

/ P'y{x)dx 
where C > is some constant independent of R. 

Remark 10. Again the same result holds true when the purely attractive Newtonian 
interaction is replaced by a potential W satisfying the assumptions in Remark\^ 

Remark 11. It is an interesting open question whether the HFB-minimizers (7, a) S 
IC\ ( or less ambitiously) the density are unique for given Tr(7) = A, modulo sym- 
metries of the HFB energy functional. 



The decay estimate for p-yix) and a slightly refined version (see Lemma 15.31 
below) will play an important role in the proof of Theorem [TJ Let us note that 
if 7 were a finite-rank operator, the self-consistent equation (written in the same 
form as in [26j ) would show immediately that the associated density p-yix) is indeed 
exponentially decaying. However, we believe that 7 is in general not a finite-rank 
operator for any q ^ 1. We shall justify this claim for the physically relevant case 
q — 2 below, provided we assume that a 7^ holds. 

3.2.2. The spin-1/2 case (q — 2). Let us now consider the specific case when the 
number of internal spin degrees of freedom is q — 2, which corresponds to the 
physically relevant case of spin-1/2 fermions such as neutrons. In this case, it is 
shown by Bach-Frohlich-Jonsson in [T], based on a concavity result in Lieb |28) . 
that the following equality holds: 



/(A) = ^(A) for q = 2 spin degrees of freedom. 



Here 

(3.6) J(A) =inf{ 2 J-(r) : t e B{l'^{R^;C)), 
and the reduced (no-spin) energy is defined as 



T, Tr(T) = A/2}, 



(3.7) TiT):=TT{TT)-KDipr,Pr) 



-Ex(t) 



y/T{l - T){x,y) 



2 ' ' 2 77r3xr3 |a; - y\ 

Furthermore, all the minimizers for £(7, a) and F{t), respectively, satisfy 
'1 0\ , , r „ / 1 



dx dy. 



(3.8) 



7 ^ 



1 



and 



±Vr(l-r) 



-1 



We note that J-{t) has a form similar to the Miiller functional which was studied 
in |14] . Indeed, a straightforward adaptation of an argument given in [14 combined 
with p. 81) leads to the following statement (its proof is outlined in Section [S]) . 

Theorem 3. (Infinite Rank of Minimizers if a 7^ 0). Assume that (7,0) is 
a minimizers for /(A) with < k < A/tt and A > 0. Then, if a ^ holds, the 
operators 7 and a both have infinite rank. 

Remark 12. Since it is expected that a ^ holds in general, this result justifies 
the variational arguments developed below to prove a decay estimate for p-y. 
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3.3. Minimizers for Simplified HF Model. As we have mentioned above, a 
rough approximation consists in neglecting the pairing term, i.e., by restricting to 
density matrices having a = 0. A result similar to Theorem [I] holds true for an 
adequately defined critical mass A"^°(k) and without the assumption that k < A/tt. 
Minimizers then satisfy the following equation 7 = X(-oo.^) (-^7) + The proof is 
indeed much simpler, as can be seen from the proof of Theorem [T] 

In the Hartree-Fock (HF) model, one restricts the search of a minimum to the 
variational set to density matrices 7 that are projectors of rank N, i.e. "/ — 
Sil=i \'Pi){'Pi\7 with some L'^-orthonormal functions {f i}fLi- For ^-ny integer A'^ ^ 1, 
it is convenient to define the set 

(3.9) Vn := {(1,0) e IC : 7' = 7}, 
and accordingly the HF minimization is given by 

(3.10) |jHF(jv) :=inf{f(7,0) : (7,0) eT'w}. 



In analogy to the HFB model discussed above, we denote by N^^{k) the largest 
integer for which I^^{N) > —00 holds. We now have the following result about 
the HF minimization problem. 

Theorem 4. (Existence of Minimizers for HF-Model). Let m > and 

K > 0. First, there is no minimizer for /^^(l). Then for every integer N such that 
1 < N < N^^{k), the following holds: 

(i) Every minimizing sequence {(7n,0)}„gN for {N) is relatively compact in X 
up to translations. In particular, there is a minimizer (7,0) S Vn for I^^{N). 

(ii) The following binding inequality holds, for all integers 1 ^ k ^ N — 1: 

(Hi) Any minimizer (7, 0) £ Vn for I^^{N) can be written as ^ = X^^^^i Wn) {fN\, 
where the functions (pn G H^^'^^E?) satisfy 

(3.11) H^Lpn = A„(/?„, 

with Ai ^ . . • ^ Ajv < being N negative eigenvalues of the the mean-field operator 
defined in ()3.3p . 



Since the HF model introduced above is an A^-body theory, one can prove, by 
using the geometrical methods of [T31 HHl 133] , that every minimizing sequence is 
compact (up to translation) if and only if the quantized binding inequality in (n) 
holds true. This was done for the first time by Friesecke in [15] in the context 
of the usual atomic Hartree-Fock and MCSCF theories. The proof of Theorem |4] 
relies on Friesecke's method, together with some ideas of the proof of Theorem [1] 
In some sense it is much easier than for Theorem [TJ since it is possible to prove 
the quantized binding inequalities by induction on A^. We briefly explain this in 
Section El below. 

3.4. Time-Dependent Theory. We now make a brief digression into the related 
time- dependent HFB theory. As it is well-known from nonlinear Schrodinger equa- 
tions, a detailed understanding of the variational calculus for the time-independent 
theory turns out to be of great use when addressing dynamical questions; e.g., the 
stability of solitary waves and blowup analysis. Here, however, we shall mostly limit 
ourselves to the dynamical (orbital) stability of solutions generated by minimizers 
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(which essentiahy fohows from Theorem [T] as a by-product). As a concluding out- 
look, we briefly comment on the finite-time blowup for the time-dependent HFB 
equations, corresponding to the onset of the dynamical collapse of a neutron star. 

For the reader's orientation, we recall that the time-dependent HFB equations 
generated by the functional £(7, a) can be written in commutator form as 

(3.12) idtT^[Fr,T]. 

Here and as usual, F = F(7,a) denotes the admissible 1-pdm defined in (j3.1l) and 
Fr is the HFB mean-field operator introduced in (13. 2[) . Clearly, the equation (I3.12p 
is a nonlinear evolution equation. With regard to the wellposedness of its initial- 
value problem, we note that a straightforward adaptation of [71 IHl IH HZl HOI US] 
yields the following results. 

Local Wellposedness. For each initial datum (70, ao) € /C and T > sufficiently 
smah, there exists a unique solution (7,0) € C"([0, T); A") n C^([0, T); A") solving 
(|3.12l) such that r(0) = F(7o,ao). Moreover, we have conservation of energy and 
expected number of particles, i. e., 

VO^i<T, f(7(i),a(t)) =£(7o,ao) and Tr(7(t)) = Tr(7o). 

Global Wellposedness. If the initial datum (70, ag) € K, satisfies the smallness 
condition Tr(7o) < A"™(k), then we have the a priori bound 

sup \\{j{t),a{t))\\x i^C 

0^t<T 

and the solution F(t) — r{j{t),a{t)) extends to all times ^ < < 00. 

Next, we consider the behavior of minimizers for the HFB energy £"(7,0) with 
respect to the evolution equation (I3.12p . To this end, let us assume that the initial 
condition (70, ao) is a minimizer for /(A) as given by Theorem[2l With this choice of 
initial conditions and using equation p.5|) above, an elementary calculation shows 
that the corresponding solution r(t) = r{'y{t),a{t)) of p.l2|) is given by 

(3.13) 7(t) = 70 and a(t) = e-^f^'^ao, 
or equivalently in a 2 x 2 matrix form 

(3.14) r(t) = e-'"'^^* Fo e*'^^* 

with /i < taken from Theorem [5] and N defined in p.4p . Hence HFB minimizers 
give rise to stationary solutions of (|3.12l) . as one naturally expects. Moreover, by 
adapting a well-known general argument in [5], the relative compactness result of 
Theorem [1] and the conservation laws imply orbital stability of minimizers. The 
precise statement can be formulated as follows. Let us define the set of minimizers 

Mx ^ {h,a)elCx a) =I{X)}, 

and introduce the distance function on /C 

distA^;,(7o,ao) = inf || (70, ao) - (7, a)|U- 

Theorem 5. (Stability of Minimizers under HFB Time Evolution). Un- 
der the assumptions of Theorem]^ the set of minimizers Mx for /(A) with < 
A < A^^^(k) is orbitally stable in the following sense: for every e > 0, there is 
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6 > such that if the initial condition satisfies dist»A (70i Q^o) < ^, then the corre- 
sponding solution r(t) = {'-f{t),a{t)) of 13.12\} exists for all times t ^ and obeys 
suPt^odistAiA(7(i),a(*)) < £• 

Finally, we briefly comment on finite-time blowup for the time-dependent HFB 
equation. Indeed, recent work by Hainzl and Schlein [21j (which extends previous 
results in [T71 [TB]) addresses blowup for the Hartree-Fock (HF) approximation. 
Since, in particular, the HFB time evolution (|3.12p preserves HF states, we can 
immediately infer from |21[ [171 the following result. 

Finite-Time Blowup for HF States. Let the initial datum (70, 0) G Vn be 

a Hartree-Fock state; see p.9p . Furthermore, assume that 70 £ C(^{M.'^ x M"^) 
and that we have radially symmetry in the sense that jQ{x,y) = "fo{Rx, Ry) for 
all R G SO{3). Then, the energy condition £(70,0) < ~mN implies that the 
corresponding solution r{t) — r(7(i),0) of (|3.12l) blows up in finite time, i.e., we 
have limt^T Tr(V— A7(t)) = -t-00, for some < T < 00. Furthermore, the solution 
exhibits minimal mass concentration at the origin in the sense that, for every i? > 0, 



where N^^{k) > is the same constant as in Theorem |4] above. 

Remark 13. In order for £ (^jq , 0) < —mN to hold, we must have large initial data 
in the sense that Tr(7o) — N > N^^(k) holds (which implies that I^^{N) = —00 
for such large N ); cf. Section \3.3\ for the definition of N^^ . 

Remark 14. It is of interest to extend the above blowup result to HFB states 
with non-vanishing pairing a ^ 0; or, even more ambitious, to relax the radiality 
assumptions on the initial data. 



In this section we collect some preliminary results that are needed to set up the 
variational calculus for the proof of the main theorems. 

4.1. Weak-* Topology on X. Recall that we assume that g ^ 1 (describing the 
spin degrees of freedom) is a fixed integer and that (3p with \ ^ p < 00 denote 
the Schatten spaces of bounded operator acting on L'^(IR'^; C). It is well-known 
that ©p are reflexive when 1 < p < 00, the dual being ©p/ where p^^ -\- p' — 1. 
Furthermore, it is a classical fact that 61 is the dual of the space of compact 
operators ©00; see [41j . This induces naturally a weak-* topology on X for which 
the unit ball is compact, by the Banach-Alaoglu theorem. More precisely, we say 
that (7,1, q;„) ^ (7,0!) weakly-* in X if 



for all K e ©00 and all K2 £ ©2- Note that the convex set of fermionic density 
matrices /C C A" is both closed in the strong topology of X and in the weak-* 
topology. However, the set IC\ is not closed in the weak-* topology. That is, 
we can have (7n,Q;„) ^ (7,0;) weakly-* in X with 7„ e IC\ such that Tr(7) < 




4. Preliminaries 



(4.2) 



(4.1) Tr((l 



A)i/V(l - Ay/^K) ^ Tr((l - A)i/S(l - A)^^^K), 
Tr((l - A)i/4a„X2) ^ Tr((l - A)'^^aK2), 
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liminf Tr(7„). Put differently, the continuous linear functional (7, a) 1—^ Tr(7) on 
X is not weakly-* continuous. 

Throughout this paper, we shall make use of the following important fact that 
helps conclude strong convergence a posteriori: Suppose (7„, a„) ^ (7, a) weakly-* 
in X and moreover assume that 

(4.3) Tr((l - A)i/V(l - A)i/4) ^ Tr((l - A)i/S(l - A)'/^), 
and 

(4.4) Tr((l - A)i/4a>„(l - A)!/*) ^ Tr((l - A)i/4a*a(l - A)^/^). 
Then we have (7„, a„) — > (7, a) strongly in X. This follows from [HI Thm A. 6]. 

4.2. Some Basic Properties of HFB Energy. We have already seen that the 
energy £ is well-defined on JC. Next, we collect some basic facts. 

Lemma 4.1. (Basic Properties of /(A)). Let m > and < k < A/tt. Then 
the following holds. 

(i) For alio < X < A^™(k), we have -00 < /(A) < and, for all < A' < A, 

/(A) ^/(A-A') + /(A'). 

Moreover, if equality holds for some A and X' , then there exists a minimizing se- 
quence for /(A) that is not relatively compact in X up to translations, 
(a) The function X i~> /(A) is continuous and decreasing. 
(Hi) The functional £ is not weakly-* lower semicontinuous on K,\. 

Proof. The proof of the subadditivity inequality in Part (i) follows from classical 
arguments: we take two states (7,0?) G ICx-y and (7', a') S ICx> that are almost 
minimizers of /(A — A') and I{X'), respectively. By a density argument, we can 
assume that (7, a) and (7', a') are both finite-rank operators with compact support. 
Next, we define the sequence of states 

(4.5) (7„,3„) = (7,0;) -I- T*g(7',Q!')r„e, with n = 1, 2, 3, . . . 

where e G K'^ is some fixed unit vector and T„e/ = /(■ — ne). Note that (7, a) and 
T*g(7', T')Tne commutc for n large enough, and we easily verify that (7„, 5;„) € 1C\ 
when n is sufficiently large. By taking the limit n 00, it is straightforward to 
conclude the subadditivity inequality. If moreover equality holds for some A', then 
(7,1, 5„) furnishes a minimizing sequence that fails to be relatively compact in X 
up to translations. 

To see that /(A) < holds, we choose a fixed state (7, 0) € K,\ with 7 smooth 
and finite-rank and such that D[p^,p~f) — Ex(7) > holds. (The later condition is 
indeed equivalent to saying that rank (7) ^ 2.) Next, we consider the energy of the 
rescaled state 75 := UslUI where iUsf){x) = 5"^^'^ f{x/S). Note that (75,0) € K-x 
for all (5 > 0, as one easily verifies. Using now the operator inequality T ^ —-^A, 
we deduce that 

£(75,0) lTr(-^A7) - ^{D{p^,p^) - Ex(7)) < 0, 

provided that (5 > is chosen sufficiently large. This establishes Part (i). 

To show Part {ii), we note that the strict monotonicity of /(•) is an obvious 
consequence of the subadditivity inequality and the strict negativity of /(•). The 
continuity of /(•) follows readily from using trial states. This proves Part (ii). 
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To prove Part [iii), we fix some < < A and consider a state (7^,0) E 
K.\i with compact support (meaning that its density p^i has a compact support). 
Furthermore, let (7^, 0) G JC\_\i be with compact support and such that £(7^, 0) < 
0, which is possible because of /a (A — A^) < 0. Define then 7„ := 7^ + T„e7^T*g 
with r„e as above. For n large enough, the terms in 7„ have disjoint supports and 
we obtain that (7^, 0) G K,\. Next, a computation shows that 

lim £:(7„,0) = f(7^0)+£(7^0) <f(7\0). 

n— >-oo 

Since 7„ ^ 7^ weakly-* in X , this strict inequality shows that £(7, a) is not weakly- 
* lower semicontinuous on X . □ 

Lemma 4.2. (Coercivity of £). Let m > and ^ k < A/tt. For ^ A < 
A^^^(k), the energy £ is coercive on JC\, i. e., we have 

£(7,a)— ?>oo as ,a)\\x ^ 00. 

In particular, all minimizing sequences for £ on JC\ with ^ A < A^^^(k) are 
hounded in \\ ■ \\x-norm. 

Proof. For any (7, a) G K.\ and 1 > e > 0, we have 

(4.6) f (7, a) ^ eTr(T7) + (1 - e)/,/(i_,)(A), 

where /K/(i_e)(A) is the ground state energy of £ on K,\ with k replaced by k/(1 — e). 
By continuity of A^^^(«;) and since ^ A < A^^^(k), we can choose e > small 
such that A < ^^^{jz^) and hence /K/(i_e)(A) > —00. Using that a*a < 7 (see 
(|2.15p ) and the strict positivity of m > 0, we see that ||(7,a)||Ar — >■ 00 implies that 
Tr(T7) 00. This fact, by (|4.6p . implies that f (7,a) -> 00 as weh. □ 

4.3. Study of an Auxiliary Functional. The HFB energy contains two non- 
convex terms: the direct and pairing term. It turns out to be expedient to introduce 
the following energy functional 

(4.7) c;(7,a) := T,{K^) ~^(( ^^^^^dxdy, 

2JJr3xr3 \x-y\ 

for (7,0:) G /C, where we set 



(4.8) K V-A + m2. 

Recall estimate p.l6p which tells us that the pairing term is controlled by the 
kinetic energy. Hence we have the lower bound 

(4.9) e(7,a)^0, 

when ^ K ^ 4/7r. We start by stating a result similar to [121 Theorem 1]. 

Proposition 4.1. Let m > and ^ k < 4/7r. Then Q is weakly-* lower semicon- 
tinuous on IC. This means */{(7n,Q^n)} C /C satisfies (7n,ctn) ^ (7; ct) weakly-'^ in 
X, then 

(4.10) liminf0(7„,a„) ^ 0(7,a). 
Additionally, equality holds if and only if the convergence is strong in A'. 
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Proof. Consider a sequence in JC such that (7,1, a„) (7, Q^) weakly-* in X. By 
Lemma [2. 11 y^p^ is bounded in H^^^{W^). Thanks to a Rehich-type theorem for 
H^^'^{MP) and extracting a subsequence if necessary, we can assume that p^^ 
weakly in LP{R^) for 1 < p 3/2 and strongly in Lfo^(M^) for 1 < p < 3/2."As for 
the sequence a„(a;, y), we notice that 

Tr((l - A)i/Va:(l - A)^/') = ^((^1 - A, + yT^)a„, a„)^,(j,3xE3), 

Thus a„(x, y) is a bounded sequence in x K.-^), and we may assume that 

an{x,y) converges strongly in L^qj,(M'^ x M'^), by a RcUich-type theorem again. 

Next, we recall the definition of the cutoff functions XR and C/j in (|l-5p . We then 
find the estimate 

\an{x,y)\^ 

-ax ay 



3xR3 \x-y\ 



C-R(a;)^|a„(a;,y)|2 ff XRix?\anix,y)\'^ ^ 

-axay+ // j j ax ay 



\x-y\ yyR3xR3 \x-y\ 

Xfl(a;)^X3fl(y)^|a„(a;,y)|2 ||a„|| 



2 



2 JjR3xR3 |a;-y| 

In the last line we have used the Hardy-Kato inequality p.l3p and the inequality 
Cfla„a*C_R ^ Cr^tiCb. thanks to (|2.15p . Because (jmCtn) is bounded in X, the last 
term in the above estimate is 0{R~^). 

Now we localize the kinetic energy. Lemma lA.ll tells us that the nonlocal operator 
K satisfies 

I 1"°° I 1 
XrKxr + C,rKCr ^K + -j^ (|Vxi?P + iVCflP) j:^^sds 

(4.11) r^K + C/R^, 

using ||Vx_r||^oo + ||VCi?||^oo ^ C/R'^. Since Tr(7„) is uniformly bounded, this gives 
us Tr(if7„) > TYiKxRlnXR) + Tr(ifCfl7«Cfl) " C/R^. Thus we obtain 

(4.12) g{-in,an) > ^T^{K{xRlnXR " XflanXsfl^nXfl)) + Tr(Xxfla«X3fi"«Xfl) 

XR{x)^X3R{yf\an{x,y)\^, , /.xrr/r^. .^ ^ 



K 



r, I I I I dxdy+{l- K7r/4)Tr(i^CK7«Cfl) - ^• 

2JJr3xr3 \x-y\ R 

Since a„ converges strongly to a in Lj^Qj,(M'^ x R^), we have that XRC^nXsR^^nXR ^ 
Xi?Q^X3KQ!*X-R strongly in 61. Hence we deduce that 

XR.lnXR - XRanxlR.a*nXR ^ Xi?.7Xfl - XflaX3fl"*Xfl 
weakly in 61. Notice also that, by p.lSp . 

XROinXlRO^nXR < XRC^n^nXR < XRIuXR 

Therefore we may use Fatou's Lemma to infer that 

Iminf Tr(ii:(xfl7«Xi?. - Xi?.a«X3fl"«Xi?)) ^ Tr(is:(xfl7Xi?^ - XRaxlnQ* Xr)) ■ 

Next, we introduce the following sequence of functions on M.-^ xM'^ given by a^'{x, y) :- 
XR{x)X3Riy)'^n{x,y), which converges weakly in H^^^{R^ x R"^) to a^{x,y) :— 
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XRix)x3R{y)a{x,y). We can write 

TT[KxRanX3R(^nXR) - 1^ , , dx dy 

^ J JW-y.M? F ~ y\ 



L2(R3xR3) 



Using that k < 4/7r, the Hardy-Kato inequality yields 

- oi I >(1 - «'^/4)^x > (1 - '«7r/4)m > 0. 
2\x-y\ 

Thus we may use Fatou's Lemma one more time to obtain 

,4.13, ,i„.„r((i..-^)„„VJ,^,^_^_^ 

R „.R 



L2(R3xR3) 



In summary, we have derived the following inequality 

r ■ rrf \^t-^&' a f f XR{x)^X3R{y)^\a{x,y)\^ ^ C 
hmmf g(7„,a„) > Tr(Xxfl7Xfl) - 7; l I \ \ dx dy - — . 



n— f oo 



2 J jR3xR3 N - 2/1 

By taking the limit i?, — > cxd, we obtain the desired result 

liminf C/(7„,a„) ^ ^/(t,")- 

n— >oo 

It remains to show the claim about strong convergence and equality. First, we 
note that if equality holds in dUTOl), then by (|4?T2l) . ^(^\^ and k < 4/7r, 

(4.14) lim limsupTr(XCi?.7„Cfl) = 1™ limsupTr(i4:Ci?,a„a*Ci?) = 0. 

On the other hand, the lower bound (14.13^ implies that in case of equality, for any 
fixed R, we have — strongly in H^/'^{R^ x M'^). We easily deduce from (|4.14p 
that strongly in H^^'^iM? x R^). By this strong convergence, we can now 

pass to the limit in the pairing term and we conclude Ymin^oo Tr(ii'7„) = T!v{K^), 
which itself eventually implies that (7„,a„) — > (7, a) strongly in X. The proof of 
Proposition 14. II is now complete. □ 



A simple but important consequence of Proposition 14 . 1 1 is the following fact. 

Corollary 4.1. (Conservation of Mass implies Compactness). Let m > Q 

and < K < A/tt. Consider a minimizing sequence {(7^, an)}neN for /(A) with A < 
A^^^(k) and such that (7„,a„) ^ (7,0;) weakly-* in X. Then (7n,Q„) — > (7,0;) 
strongly in X if and only ii/Tr(7) = A. 

Proof. Clearly, if 7„ — 7 strongly in X , then Tr(7) = A. To show that Tr(7) — 
A is also sufficient for strong convergence, we argue as follows. Recall that any 
minimizing sequence (7„,a„) with A < A^^^(k) is bounded in X, by Lemma 14.21 
Then Lemma [2TT] implies that is bounded in H^/'^{M?). Thus we may assume 

that weakly in H^/'^{M?) and p^^ — >■ pointwise almost everywhere. 



MINIMIZERS FOR HFB THEORY OF NEUTRON STARS 



19 



Using this pointwise convergence and that p^^ < C, we deduce from the Brezis- 
Lieb refinement of Fatou's lemma (see, e.g., |29l Theorem 1.9]) that 

p^^{x) dx — / p^{x)dx+ / \p^{x) — p^^{x)\ dx + o{l) , 

where o(l) — >■ as n — oo. Therefore Tr(7) = J pj = X implies that p-y^ — > p-y 
strongly in L^(R'^). Furthermore, by interpolation and the fact that |1p7„|1lj' < C 
for l^p^3/2, we deduce that 

(4.15) pj^ ^ p-y strongly in LP{m.^) for 1 «C p < 3/2. 

Thus the Hardy-Littlewood-Sobolev inequality shows that D{py^ , p^^) — > D{py, p^). 
The exchange term is found to be weakly-* lower semicontinuous by Fatou's lemma: 

(4.16) liminf Ex(7„) > Ex(7). 

By Proposition 14. 1[ we have that the kinetic energy plus the pairing term form a 
wise-* functional on 1C\. Hence, we conclude 

(4.17) liminf£(7„) ^f(7) ^/(A). 

n— f oo 

As by assumption £'(7n) this proves a posteriori that 7 is a minimizer for 

/(A) and that 

(4.18) lim g(7„,a„) = g(7,a). 

By Proposition 14. II again, this shows that (7,1, a„) — ?> (7, a) strongly in X. □ 

Next, we study the minimization of the functional Q in more detail. This will be 
useful for the study of a minimizing sequence for /(A) when ruling out vanishing; 
see Section [71 We start by introducing the following variational problem: 



(4.19) G(A) := inf{g(7, a) : (7, a) S /Ca}. 



Indeed, we can derive the following formula for this infimum. 

Proposition 4.2. (Value of G(A) and upper bound on /(A)). Let m > and 

< K < A/tt. Then for ^ A ^ A^^^(k), we have the equality 

(4.20) G(A) = /3A, 

where m{l — ktt/A) < /3 < m is given by 

inf {ib,H^) ifq^l, 

V'eiodd(R')JI'ML2=i 

inf {ip,H'iP) tfq^2. 

^GL2(K='),||V'lli2=l 



Here the integer q ^ 1 denotes the spin degrees of freedom, and H = V— A -|- 
is defined in the sense of quadratic forms on L'^{S?) with form domain H^^^ 

Moreover, we have the strict inequality, for all < X < A^^^(k), 

(4.21) /(A) < G(A) - mA. 

The proof of Proposition 14.21 is given in Appendix O It is partly inspired of the 
discussion of the Miiller energy functional in [14]. 



Remark 15. Adapting the proof of Proposition IJT^ given in Avvendix[^ one can 
show that G(A) = /(A) ~ —00 if k > A/n. Indeed, in this case the one-body operator 
H is not bounded below anymore; see [22j . 
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5. Proof of Theorem [2] 

In this section we prove Theorem [2] which estabhshes key properties of any 
minimizer (7, a) of the HFB energy functional £. In fact, we shah need some of these 
properties (more precisely, a certain decay estimate) in the proof by contradiction 
of Theorem [1] that establishes the existence of minimizers; see Section [7] below. 

5.1. Mean-Field Equations and Bounds on the Lagrange Multiplier fj,. 

Let (7, a) G /C a minimizer for /(A). To derive the equation satisfied by (7, a), 
we can follow a general argument given in [5]. That is, by writing that £((1 — 
t)(7,a) + t(7',a')) > £(7,0;) for all {j',a') S IC\, we see that (7,0;) is a solution 
of the following linearized variational problem 

(5.1) K{X) = i inf Trpo {Fr - P^)} , 

2 (7',a')e/C " 
Tr(7') = A 

where Trpo (•) denotes the P^-trace as defined in We use the definition of Fr 
and r introduced in p.2p and (|3.ip , respectively, and we employ the notation 

Note that 

(5.3) Tr(7') = ^Trpo {iV(r'-Pi')}, 

and that the function A 1— > K{X) defined in (|5.ip is easily seen to be convex. Hence 
K{\) satisfies 

(5.4) K{\') ^ K{\) + ^(A' - A), for all A' > 0, 

whenever ^ G [_fi';^(A), _ft';^(A)], where K'_{X) and K'^{X) denote the left and right 
derivatives of K{X), respectively. (Note that K'j.{X) indeed exist, since A 1— K{X) 
is convex.) Following the proof of Lemma 14. 1[ we find that K{X) is non-increasing, 
and thus ^ ^ holds. Moreover, by using (|5.4p . we deduce that (7, a) is also a 
minimizer for the following problem without constraint on the trace: 

(5.5) ^inf Trpo {(Fr - fiN) (F' - F^) } , 

where A'' is defined in p. 41) . This observation allows us to deduce (like in [TH]) that 
F = F(7, a) solves the equation 

(5.6) T = X(-oofi){Fr-t^N)+D, 

where Z? is a finite rank operator of the same matrix form as F and satisfies ran(Z?) C 
ker(Fr — fJ-N). As an important next step in our analysis, we now prove an upper 
bound on the chemical potential /i . 

Lemma 5.1. (Upper bound on /i). Let F = F(7,q:) be a minimizer for /(A) 
and let ^ be as above. Then we always have that < holds. If moreover a 7^ 0, 
then fj, < P — m holds, where /3 < m is given by Proposition \4-.S\ above. 

Proof. The proof in the case a = (Hartree-Fock case) is standard and deserves 
no further discussion. Hence we will assume that a ^ from now on. We start by 
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expanding the energy around the minimizer F up to second order. We make use of 
Bogohubov transformations of the form e"^ with 

where h G &i and k E &2- From |2] we recah that, for every e g R, we have that 
Fg — e^'^Fe^"^ is an HFB state. Expanding up to second order, we find 

Fe = F + eFi+e2r2+0(e3), 

where 

ri=i[fl",F], T2 = HTH -^H^r -^TH^. 
Next, we take h = and k = —ia with a = —a^ E ©2 to be chosen below, that is: 



H 



—ia 
ia* 



An elementary calculation yields 

"^ I /^7i "1 \ r ^ ^ \ + {^'^ "2 



with 



— aa + aa , ai = —7a — aj, 
1.1. .1 



72 = —070 — —aa 7 — 2'^'^'^ I '^2 = —aa a — —aa a — ^'^^ 
Calculating the energy and using that [Fp — /lA^, F] — 0, we find 
£{T,) = S{T) + eMTr(7i) + e'Qia, a) + O(e^), 

where 



Q(a, a) — Tr{T{a*a + 72)) — kD{p~^, pa*a+-t2) + k3?A(7, a*a + 72) — k3?X(q;, 02) 

- 1^(^71 + 1^(71,71) - |A:(a + ai,a + ai). 
Here we use the convenient notation 



(5.7) Xia,h)-ll ^^^p^^dxd,. 

JJwxM? \x-y\ 

We now claim the following fact. 

Lemma 5.2. There exists a — —a^ with TT{Ka*a) < 00 such that Tr(7i) < and 
Q{a,a) < (/3-m)Tr(72). 

Assuming Lemma 15.21 for the moment, we have Tr(7e) < Tr(7) = A for e > 
small. Furthermore, with the help of Proposition we deduce 

£-(F,) ^ /(Tr(7,)) ^ /(A) - /(A - Tr(7,)) ^ /(A) - (/? - m)(A - Tr(7,)) 

= /(A) + (/3 - m)(eTr(7i) + e^Tr{j2) + ©(e^))- 

Collecting these estimates, we obtain 

(/3 - TO - M)Tr(7i) e{Qia, a) - (/3 - m)Tr(72)) + O(e'), 

whence the desired estimate on /i follows, by taking e > sufficiently small. This 
concludes the proof of Lemma |5. 11 provided that Lemma 15.21 holds, which we will 
show next. □ 
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Proof of Lemma \5.2l We begin with a preliminary observation; Assume that we 
have constructed an a with = —a and Tr{Ka*a) < oo, such that Q{a,a) — 
(/3 — m)Tr(72) < 0. If Tr(7i) > 0, it suffices to replace a by —a to conclude the 
proof. If Tr(7i) = 0, then we can replace a by a' = a — ija. We then have (with 
an obvious notation) Tr{j[) = -2riTT{aa*) < and Q{a',a') - {(3 - TO)Tr(72) = 
Q(a, a) — (/? — m)Tr(72) + 0{r]). Hence we can make the last term < 0, by choosing 
r] > sufficiently small. As a conclusion, we only have to construct a such that 
Q{a, a) < {(3 — m)Tr(72) without considering the constraint Tr(7i) < 0. 

To construct a, we use trial states in the spirit of Appendix [Cl below. Here, let 
us for simplicity assume that we have q — 1 spin degrees of freedom (since the proof 
for g ^ 2 is analogous). That is, we choose 

aR,L{x, y) = Xl{R{x - 3Lv))f{Rx - Ry)xL{R{y - SLv)) 

where R G SO{3) is a rotation and ?; G M"^ is a fixed normalized vector; and 
/ S iJ^/^(K^) is the normalized ground state solution to (K — k/2|x|)/ = /?/, with 
K — k/2\x\ acting on ^^^^^(M^). Moreover, for L > given, we define xl{x) ~ 
L^^/'^xix) with X S (^^(R"^) being a radial function with support in the unit ball 
{|a;| ^ 1}. We will use the following fact: 

C 



(5.8) 



L3/2' 



where || • || denotes the operator norm on L^. Assuming this bound, it is straight- 
forward to derive that 

C 



71 



2,3/2' 



Tr(if72 



72 



6i 



1?' 



Using these bounds and the fact that \X{f,g)\ ^ C||i4r^/2/|| 
deduce that 



we 



Q{aR,L, aR,L) ~{f3~ m)Tr(72) < Ti{{K - /3)a,j,ia^^i) 



-j^X{aRx,aRx) 



The first term of the right-hand side is invariant under translation and rotations 
(and thus independent of i? G SO{Z)). The estimates in Appendix [Cl give 

Ty{[K ~ f3)aR,La*n^i) - ^X(afl,L, afl,^) = O {L-^) , 

which implies 

Q{aR^L,aR,L) - (/? - m)Tr(72) s=: -^D{p^, p^^.a-^J + O (l-^"') . 
Now we average over R E SO{3) and apply Newton's Theorem to obtain 

dR D[p^,pai^^^a% 



SO(3) 



> Tr(ai,Lat.L)/|:.KLP7 e 



for L > large enough and some constant e > 0. Hence the conclusion of Lemma 
5.21 follows for a suitable rotation R — R{L). 

It remains to prove the key estimate (15.81) . First, by rotational and translation 
symmetry, we can assume that i? = 1 £ 50(3) and w = G in the definition 
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of aR^L above. Thus it suffices to consider aL{x,y) — XL{x)f(x — y)xL{y) and to 
show that ||iir^/^aL|| ^ C/L^^^. To prove this, we note that ql acts hke a^tp = 
XlU * (XiV")): for V' € L^(R^). Next, we observe that 

(5.9) K'^'aLi^ = xlU * (XL0) + XlUK'^'J) * (xl^)) 

Since f{x) decays exponentially (which, e.g., follows from the method in |10j). we 
conclude 

(5.10) \\K'/\f\\L^^C. 

Indeed, the first bound follows immediately from the exponential decay of f{x). To 
derive the second bound, we rewrite the equation satisfied by / to find that 

Note that the integral kernel of K^^"^ = + m'^)^^^'^ with m > belongs to 
i^(M^); see, e.g., |42l p. 132]. Since + PJWl^ < oo (by exponential decay 

of /), the second bound in (|5.10p follows. Going back to (|5.9p . we now obtain that 

\\K'/W^L^ ^ ll./llLi||[i^'/',XL]|l||XL|lL~||0|U^ + WK'^'fhAlXLWl^mL^ 

«;c(L-'^/2^i-3/2)||0IU2. 

Hereweusedthat ||[i^i/2,XL]|| ^C\\^Xl\\l^ <:CL-^-^/^ and ||xl|| = ||xl||l- 
CL^^I^. This implies the bound (j5.8p and completes the proof of Lemma [521 D 

5.2. Decay estimate. The next step consists in proving the following decay result 
for the density function of minimizers. 

Lemma 5.3. For all R > sufficiently large, we have 

C 



(5.11) TriTCR-fCR) + J CUx)P',{x)dx 

where (^r is the cutoff function defined in (|1.5p and C > denotes some constant 
independent of R. 

Proof. When a = 0, we have /i < by Lemma [??T1 hence 7 — X(-oo,o) (^7 ^ A*) + ^ 
is finite rank. Exponential decay is obtained by following the method in |10j . 

We now assume that a 7^ 0. By Lemma TS. 11 we have fx < [3 — m. Writing that 
{Fr — fJ.N)r ^ and extracting the first diagonal term yields 

(5.12) A:^ (H-y ^ fi)-f - nX^a* i^O, 

where Xa denotes the operator with kernel Xa{x, y) ~ a(a;, y)l\x — y\. 

On the other hand, we consider [Fr — /iA^, F] = and extract the upper right 
corner. This yields the following two-body equation for the pairing wavefunction 
a = a{x, y) G H^/'^{R^ x R^). 



(5.13) yiH-r)x + {H^)y - - 2^j a = -K{-f^ + 7y)^a, 

Here means that A as an operator on L^(R^) acts on the x-variable in a{x,y). 
Likewise Ay means that A acts on the y-variable in a{x,y). Our goal is now to 
combine (|5.12p and (|5.13p in a useful way to extract decay information. 
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Step 1: Estimates from KU^ . Writing Tt{Cr{A + A*)Cr) < where A = A* 
is defined in ()5.12|) . we deduce that 

(5.14) Tr I 7 j -t^J^^ Cr{x)P-({x) dx 

[[ CUx)^-^^^^dxdy + K f V{x)CUx)Py{x)dx, 

jjk3xr3 f - y| J 

where we put V{x) = Note that we discarded the contribution of the 

exchange term in Hj, because it gives a nonnegative contribution to the left-hand 
side of the above inequality. 

Since y^p^ e H^^'^{M.^), one can see that V{x) as |a;| — > oo. Hence the term 
J V(ffP-f can be estimated by S{R) J (j^Pj where d{R) = sup^,-^^^j^\V{x)\ — > as 
i? — > oo. Moreover, the localization formula (jA.2l) from Lemma [A . 2 1 vields 

(5.15) Tr(rCfl7Cfl) - M / CRi^)pyix) dx ^ S{R) [ (:l{x)p^{x) dx 



+«// CM^) ' ; i' dxdv+-i:,[ / ^^ivCflp-^yids 7 . 

The last term can be estimated as usual by ^ CXjR}. Using that in particular 
/i < 0, we finally deduce 

(5.16) Tr(TCfl7Cfl) + / Clkx)p^{x)dx ^^+C [[ cU^)^-^^^^ dx dy , 

with C > some constant and for R > sufficiently large. We will now estimate 
the right hand side of (|5.16p by using the two-body equation (|5.13p . 

Step 2: Estimates from (I5.13p . We multiply equation (|5.13p by CR{x)CR{y) from 
the left and we project onto 

(5.17) aR{x, y) = C,R{x)a{x, y)C,R{y). 

This gives, where (•, •) is the inner product on L^(R'^ x R^), the following 

{aR,{Tx +Ty)aR) ~ K [ [ \^Ri^^y)\ dxdy -2p [ [ \aR{x,y)\'^ dx dy 

JJR3xR3 \X-y\ JJR3xR3 

= 1 + 11 + 111 + IV, 

where 

I = {aR,[T, + TyXRix)CRMa), ^^^'^''jj V{x)Cl{x)Cl{y)\a{x,y)\'' dx dy, 

III = -K{aRXRix)CRiy)iiXj)x+{X^)y)a), IV = K{aR,CR{x)CR{y){lx+Jv))Xa)- 
First, we note that II is easy to estimate by 

(5.18) \II\^S{R) f Qlix) p^(x)dx, 

JR3 

with 5{R) — > as i? oo, by using that paa* ^ p^ and F — > as — > oo. 
Next, we claim that 

(5.19) \III\ + \IV\ ^ 5{R) ( [ CUx)Py{x) dx + Tt{TCrjCr) 

\JR3 
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For ///, this can be seen as follows (we only consider the term in /// that contains 
{Xj)x, the other being treated similarly): 

(5.20) \{aRXRix)CRiy)iX^)xa)\ = \Tria*j,K'/^K-'^^Xi^^aCR)\ 

||ii:^/^Q;fl||62ll^"^''^-''^Cn7lle2l|aCj?Jle2, 

using that ||v4||eoc ^ Il^lle2 holds. Next, by the Hardy-Kato inequality, we have 
that \\K-'^/'^Xf\\e^ C\\K^^^f\\e2- Therefore we find (using that aa* ^ 7 and 
7^ ^ 7 and Cfl; ^ 1) the following bound: 

\^m\ ^ C\\K'/^an\\eM^^^CR7\\eM(:B.\\e, 



s^C^Tr(ifCflaflCfla|,Ci?)VTr(i^Cfl7'Cfl)VTr(^C/?«fl«flCfl) 

CTr(XCfl7Cfl)'/' m) Ci(a^)P7(2^) dx + Tr(TCj^7Cfl)) , 

where d{R) — ^ as i? — t- 00. The proof is the same for IV. 

Finally, we estimate / as follows (we again only consider the term with T^): 

[r^,Cfl.(2;)]aCi?(y)>l =^ \\[Tx,CR]\\\\aB.\\e2\\aCB]\e2 

C C f 

^ ^Tr(Ci?,7Cfl) CR{^)p-f{x)dx, 
n n Jr3 

where we used the commutator estimate || [Tj., /(a;)] || ^ C||V/||ioo combined with 
the fact that aa* ^ 7. 

In summary, we have proved the following estimate 

(5.21) {an, (T, + T,)afl) -k ff M^'Vf dx dy 

JJmxm \x-y\ 

\aR{x,y)\^dxdyi^S{R)(f (:l{x)p^{x) dx + TriTCRjCR.)) ■ 

Next, we turn to the left-hand side and derive a lower bound as follows. Recall 
from Appendix [C] the lower bound 

1 

'\x~y\, 

for all a G H^/^{M.^ X"^) A H^^^{M.^ ,Ci), provided that k < i/ir holds. Note that 
/? = /?('*) depends continuously on k. Using now the strict inequality 2/i < 2(/3— m), 
we deduce that 2/Lt < 2/3(k + e) — 2m for e > small, and hence 

(5.23) {aR, in + Ty)aR) ^ (k + e ~ e) [ [ M^'y^\^ dx dy 

JJr^xr^ \x-y\ 

-2fi[[ \aRix,y)\'dxdy^e If M^^^d^dy, 



(5.22) {a,[T, + Ty~Kj-—^]a) ^ 2(/3 - m)(a, a) 



\x - y\ 

for some e > sufficiently small. Therefore, we arrive at the following estimate 

(5.24) // awci(.)i.(.,.)p,^,^ 

■/■/m3xr3 \x-y\ 

^ HR) ( I Cl{x)p^{x) dx + Tr(^Cfl7Ci^)^ 
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with 5{R) — > as i? — > oo. 

Step 3: Combining (|5.16p and (|5.24p . We define the sequence {Rn}'^=i of radii 
given by i?„ =4", for n ^ 1. From (|5.16l) we conclude for all n large enough that 

(5.25) Tr(rCfl„+,7Cfl„+J + / CL^,{x)' p^{x) dx 

JR3 



2 



+ C// Cl^. W ; i' dxdy. 



(i?„+i)2 yyK3xR3 ia;-yi 



Next, we let Xi?„ = y 1 — C^^ so that C^,^ + x1e„ = 1 holds. Therefore we can split 
and estimate the pairing term in the inequality above as follows: 



// 



+ \x-y\ JJr^xr^ ^ \x-y\ 

Ck,^i-)Ckiy)n^^dxdy^:I + II. 



\x - y\ 

Since i?„ — jRn+i, we deduce from the support properties of Cr„+i (x) and Xfi„ (y) 
that ^ 

using Paa*(x) ~ /jj3 |a(x,2/)p(iy ^ Pi(x) thanks to aa* ^ 7. Further, we notice 
that Ci?„+i =^ Cfl,i and thus 

II C:lSx)(:lSv)^4^^dxdy < 5(i?„) ( I eR,Xx)p, + Tr(TCi^„7CflJ 



\x-y 

where we used inequality (I5.24p in the last step. 

Let us now define the sequence {wn}^i of nonnegative numbers given by 

(5.26) ur,= [ ClAx)p^ix)dx + Tr{TCR,,-f(:Rj- 

Collecting the previous estimates, we obtain the recursive inequality: 

C 

(5.27) Un+i ^ dnUn + -— for n ^ no. 

[Rn+lY 

where 5„ — ?> as n — )■ cxd and no 1 is sufficiently large. A simple induction 
argument then shows that u„ satisfies the following bound: 

(5.28) w„ — for n ^ 1, 

[RnY 

where _B > is some sufficiently large constant. 

It remains to extend the bound in Lemma to all i? > sufficiently large. To 
this end, let i? ^ 4 and n ^ 1 such that 4" s$ i? < 4"+^ From ([5T1| and the fact 
that ^ C4"-!; we have 

Tr(TC^7Cfl) + / eR{x)p,{x) dx^-^ + C II Cl , (x)^^^i^ dy. 

Our previous estimates show that 



Ci?,„_i(a;)-^ — ^ — 7- dxdy < 



|a;-2/| ^ 42(n-i) ^ r2 ' 



MINIMIZERS FOR HFB THEORY OF NEUTRON STARS 



27 



and the desired result follows for i? > sufficiently large. The proof of Lemma 15.31 
is now complete. □ 

Remark 16. One can bootstrap the decay estimates obtained above by commuting 
jVCflP with {K"^ + s)^^ on the right side of (IS.lSp . leading to a bound Ck/R^ 
for all k ^ 1 in (|5.11|) . This better decay estimate is, however, unnecessary for our 
existence proof, and hence we will not give any details here. 

Having established Lemma [5.31 the proof of Theorem [2] is now complete. ■ 

6. Proof of Theorem [3] 

First, we recall the equality (j3.8p and we note that, by assumption, we have 
that r(l — t) 7^ holds. To show that r must have an infinite rank, we note 
the similarity between J-"(r) and the Miiller functional studied in [T3]. Indeed, by 
following an argument in [14] . we can prove that r has infinite rank as follows. For 
the reader's convenience, we provide the details of the adaptation. 

We argue by contradiction and assume that r is finite rank. Hence we can write 

K K 

T = '^ni\ipi){ipi\, < rij ^ 1, n., = A/2 

i=l i=l 

with Tpi = ifi normalized in L^(M'^,R). Moreover, we can assume < ni < 1 
without loss of generality since a ^ by assumption. 

Lemma 6.1. We have e i7'*(M^,M) for alli = l..K and all s^O. 

Proof of Lemma \6.1i Note that tpi, ipK solve the following minimization problem 

(6.1) inf<^ ^n,(r^„¥.,)-« J J A Ldxdy 



K 



—dxdy—— 1 1 j j — dxdy, 



2JJ \x-y\ ^ 2yy \x-y\ 

((^1, ¥>k) e i^^/'(K^ M), {^,,ip,) = 6,, 

Hence they must solve the following nonlinear equation, 

.Y.k=i^k'fk{x)fk{y) 



(6.2) n,(^^2K(^^nkipk{x)^*-^^+K 



\x - y\ 



^-ni Y.k=i - nk)tpk{x)ipk{y) 



n^ \x-y\ ^ 

where the \ij are Lagrange multipliers associated with the constraints {tfi,Lpj) — 
Sij. The result is then proved using > and a simple bootstrap argument. □ 

Next, let e > be small enough and consider a real function ^ G span((pi, i — 
l..K)-^. We introduce the following test state 

K 

T, = (ni - e)\ifii){ifii\ +'^n,\ip^){ip^\ +e|'0)(V'| 

i=2 
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and compute 

(6.3) Hr.)-Hr)-^^e[[ «(^. WC^) ^ q^,)^ 



where a :— ^ r(l — t) 7^ is a (real) finite rank symmetric operator acting on 
i^(M^). Hence we will get a contradiction once we can prove the following: 

Lemma 6.2. There exists ip E span((pi, i = l..K)-^ , t/j = tp, such that 



a{x,y)ip{x)ip{y) 



dx dy = AttS^ \/ni(l - n^) f 



y\ ^ Jv? W 

Proof. Assume on the contrary that 

a{x,y)'il){x)'il){y) 



-dk > 0. 



// 



-dxdy — 



\x - y\ 

for all i/j £ span((y9i, i = l..K)-^. This means that the (real) symmetric Hilbert- 
Schmidt operator R with kernel a(x,y)\x — y\~^ vanishes when it is restricted to 
span((/?i, i ~ l..K)-^. As i? ^ 0, this implies that we must have 

(6.4) "''^^1 = 

I"" ^1 ^,J = 1 

where (aij) is a real non- negative symmetric matrix. Multipliying (j6.4p by |a; — ?/| 
and taking x = y we deduce that a{x,x) = X^ili \/'T'i(l — ni)ipi{x)'^ — for every 
X € R-^. This implies a = ^r(l — t) — which is a contradiction with ^ r. □ 

The proof of Theorem |3] is now complete. I 

7. Proof of Theorem [T] 

In this section, we present the proof of Theorem [1] By Lemma 14.11 we note 
that Part (ii) will be an immediate consequence of Part (i). Hence it suffices to 
prove Part (i), where we argue by contradiction as follows. We suppose that < 
A < A^^^(k), and we assume throughout the proof that there exists a minimizing 
sequence {(7™, Q;n)}„gN C IC\ for /(A) that is not relatively compact in X up to 
translations. We divide our proof into several subsections as follows. 

Step 1: No Vanishing. We begin by ruling out vanishing of the sequence {7n}neN- 
By vanishing, we mean that the associated sequence of densities {p-y„}neN satisfies 



(7.1) limsup sup / p^^{x) dx \ —0, 

for all R > 0. For each density matrix 7„, we define the following concentration 
function i34i 



(7.2) 



Q\{R) = sup / pj^{x)dx, 

yeK.^ J\x-y\r^R. 



Furthermore, we introduce the number 

(7.3) A^ := lim limsupQ^(i?), 
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Indeed, we shall rule out that Ai = holds. This implies that the sequence {7n} 
cannot vanish in the sense defined above. 

Lemma 7.1. (No Vanishing). We have > 0. 

To prepare the proof of Lemma FTTI we first establish the following fact, which is 
an adaptation of a similar result in [341 Lemma LI] to the fractional Sobolev space 
^i/2^r3-j ^j-^ essential ingredient will be the localization formula in Lemma [A. II 

Lemma 7.2. Let {/„} be a bounded sequence in i7^/^(M'^) and suppose {\fn\^} 
vanishes in the sense that 

limsup sup / l/n(a;)P dx \ — 0, for all R> 0. 

n-i-oo yyem^ J\x~y\^R J 

Then |1/„|Up ^ /or 2 <p< 3. 

Proof. We consider a smooth partition of unity {xfelfegz^ with J2kez^ xt = ^ ^^'^'^ 
that Xfc = 1 on Ck and suppx^ C C^, with the half-open cubes Ck = [fc, 1)'^ and 
C'f. = [fc — 1, k + where k G Z^. Furthermore, we assume that supj^g^a |Vxfc(a;)| ^ 
C. Next, we estimate as follows. 



/ \f^{x)f/'dx=J2 f Xk{xf\fn{x)f/Ux 

= E / \tci{x)fn{x)\^/^\Xk{x)f^{x)\^dx 



X/ ll-^-cJ/nl llLi(K3)||Xfe/n||^3(R3) 



(7.4) ^ C sup ||lcJ/nlllwR3) E (fn,XkKXkfn 



We have used the Sobolev-type inequality ||/|||3(r3) ^ C{f,Kf). By Lemma IAT| 
the nonlocal operator K satisfies 

(7.5) J2 ^^^^^ "^^^ir -Vho { E '^^^ ^^ds^K + C\ 

since the integral expression on the right side in the first inequality is a bounded 
self-adjoint operator due to the fact that X^fcez^ \^Xk{x)\'^ G L°°(M'^) by our choice 
of the partition {xfcjfcgz^- Therefore, we conclude 

(7-6) ^ {fn,XkKXkfn) ^ C||/„||^i/2(R3)- 

feGZ3 

Since, by assumption, we have||/„|| jLfi/2(i{3) ^ C independent of n, estimate (|7.4p 
leads to 

(7.7) / \fnix)\'^'dx ^ C sup ||lcJ/n|llf(R3) ^ 0, 

using that {|/np}neN vanishes. This shows that /„ — > in L*/'^(M'^), whence it 
converges to strongly in LP(R^) for all 2 < p < 3, by interpolation and the fact 
that {/njnsN is a bounded sequence in L^(R'^) ni^(R'^) due to Sobolev embedding. 
This completes the proof of Lemma 17.21 □ 
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We are now ready to prove Lemma l7. II 

Proof of Lemma \7.1\ By the coercivity of £ stated in Lemma 14. 2[ we know that 
our mmimizing sequence (7„,a„) is bounded in X. Hence we deduce from Lemma 
12.11 that y/]}^ is bounded in H^/^{M.^). Applying Lemma 17.21 to the sequence 
fn = ^/p^ , we obtain that strongly in LP{M.^) for 1 < p < 3/2. By the 

Hardy-Littlewood-Sobolev inequality and (|2.1ip . this implies 

(7.8) D{p^,^,p^J ^0 and Ex(7„) ^ 0. 
Therefore, we conclude 

(7.9) £{jn,an) = e(7n, ^n) - mX + o(l) ^ G(A) - mX + o(l) 

where we recall that Q and G were defined in Section [4.3l But we have f (7^, ctn) ^ 
/(A) and /(A) < G(A) — mX, by Proposition 14.21 Thus (17.91) yields a contradiction 
and we must have A^ > 0. □ 

Step 2: Dichotomy and Local Compactness. We already know that our min- 
imizing sequence {(7n,an)}nGN with Tr(7„) — X does not vanish and we have 
assumed that it is not relatively compact in X up to translations. By Corollary 
14.11 this means that we must have < A^ < A. We will now use the following 
adaptation of a the classical dichotomy result [31] to our operator setting. The 
precise statement reads as follows. 

Lemma 7.3. (Strong Local Convergence). Up to extraction of a subsequence, 
there exists a density matrix E IC with Tr(7"'^) = A""^, sequences {i?„}„gN C M+ 
with Rn 00 and {yn}neti C M'^, such that 

(7.10) y/KTlx{Rn'-)lnX{Rn'-)ry^^ ^ Vk^'Vk 
weakly-* in &i, where K = A + m^, and 

(7.11) rlx{Rn'-)lnX{Rn'-)ry^ 7' 

strongly in &i. Here Ty is unitary action representing translations in K"^ as defined 
in (II. 6p . and ^ x ^ 1 the smooth cutoff function introduced in (jl.Sp . Moreover, 
we have 
(7.12) 



lim 



Proof. We will not detail the proof which is an adaptation of ideas by Lions [31] , 
where one introduces another sequence of concentration functions 

(7.13) Ql{R)=siip p^ rj^{x)dx 

in addition to the sequence defined in (|7.2p . We refer, for instance, to 

[15] where a similar argument has been detailed. The rest of the proof follows 
well-known ideas of Lions |34| coupled with the fact that our minimizing sequence 
{7„} is bounded in X, hence (up to a subsequence) converges strongly locally in 
the trace class. The proof is the same in A". □ 
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As we have already mentioned, our model given by f(7,ck) is invariant under 
spatial translations. Thus we may assume that (in Lemma 17. 3|) the sequence of 
translations is given by 



y„ = 0, for all n ^ 1 



which we shall do for all the rest of the proof. 

Step 3: Splitting of the Energy. The next step shows that we must have 

(7.14) I{\)=I{X^) + I{\-X^) and /(A^) has a minimizer. 

To streamline our notation with Lemma 17.31 we introduce a new radial cut-off 
function ^ x ^ 1 such that x = 1 for \x\ ^ 7/5 and x = for |a;| ^ 8/5. Given 
the sequence {i?„} from Lemma [7.3[ we define the functions Xn — xi'/^n) and 
Cr„ = - xl- Likewise, we define {(7;^, )}„gN and {{'y^,al)}neN by 

ill^ai) Xn{ln,an}Xn and (7^, a^) := Cfl„ (7n, ari)Ci?„- 

Note that, since (7„,a„) G JC, one easily verifies that (7,\,a'„) e /C for i = 1,2. 
Furthermore, we have 7^ — 7^ in ©1, by Lemma [7.31 which in particular implies 
that Tr(7i) ^ J p^^ix)xlix) dx ^ A^. 
Next, we shall prove that 

(7.15) £(7„,a„) ^f(7^,«i)+f(7',«')+o(l)- 

To deal with the kinetic energy, we use Lemma I A . 1 1 and Tr(7„) = A to find that 

(7.16) Tr(r7„) ^ Tr(T7i) + Tr(T72) - CXR^'. 

The next step in the proof of (|7.15p is to separate the direct term, which can be 
done as follows: 

DiPj,, , P7„ ) = D{p^i^ , p^i ) + D{p^2 , ) + 2 J J ^'^"^J^'^yl^^^ dx dy. 
Here, we note that 

J J Xn {xfpj„ {xKn {yfpt„ jy) 



dx dy 

y\ 

Xn(xf p^Jx)C,n(.yf p^Sy) 



\x\^6B.,j5 J\y\^7B.n/S' \X — y\ 

Xnixfp^^ {x)Cn {yfp^^ (y) 



dxdy 

dx dy 

6i?.„/5^|a;K8H„/5 J\y\^9R„/5 \X — y\ 

Xn{x)'^P^„ {x)Cn{yfpj„ {y) 

6i?„/5<|2:|sS8i?„/5 J7fl:„/5^|iyK9i?„/5 \x — y\ 



Therefore, 

II 2 



Xn{x)^P^AxKniy?P^Ay) ,,^^ lOA^ i 

_ y| ^ ^ JJ^ ^ \\f^1^^\\Le/^(,B{0,9R„/5}\B{0SR„/5)} ' 

Furthermore, by (|7.12l) . we have that Pf„'^B(o,2R„)\B(Q,R,t) ^ in L^(M'^). Since 
this is also a bounded sequence in L^'(B(0, 2i?„) \ i3(0,i?„)) for all 1 ^ p ^ 3/2, 
we infer by interpolation that it also tends to in L^'(B(0, 2i?„) \ 5(0, i?„)) for 
1 < p < 3/2. Hence, by the Hardy-Littlewood-Sobolev inequality, we conclude 

D{Pt„,PiJ = D{p^i,p^i) + D{p^2,p^2) + o{l). 
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As usual, the exchange term Ex(7) is easily estimated by noticing that Ex(7„) 
Ex(7^) + Ex(7,^J. The pairing term is treated similarly as done before several times 

7.17 / / ■ ■ dxdy ^ — Tr a„a„) 

J J \x-y\ Rn 

Xn{x)'^Cn{y)'^\an{x,y)\'^ ^^^^^ 
6R„/5'^\x\^8R„/5 J7fl„/5^|yK9i?.„/5 \x — 2/| 

Note that 

(7.18) // |Q;„(a;, y)pda;dy = Tr(a„a*) < Tr(7„) = A, 

by (|2.15p . Define rjn = ri{-/Rn) for some smooth function 77 such that 77 = 1 on the 
annulus {6/5 ^ |a;| ^ 9/5} and 77 = outside the annulus {1 ^ |a;| ^ 2}. The last 
term of (|7.17p can be bounded above by 

Xnix)^Cn{y)^\an{x,y)\^ 

. . dxdy 

\x - y\ 



6i?.„/5sS|a;K8i?.„/5 "'7fl„/5^|yK9i?.„/5 



< 2^'^{K'qnan'nlalr]n)) ^ -Tr(i^7;„7„?7„) . 

Finally, the last term tends to zero by ()7.12p and thanks to the estimate || [Vk, rjn] \\ ^ 
CR^^. In summary, we have derived the estimate ()7.15p . as desired. 

Next, we notice £{-fl_,al) ^ /(Tr(7^)) and £{-f^,al) ^ /(Tr(72)). On the other 
hand, by Lemma |4. 11 

/(A)s^/(Tr(7^)) + /(Tr(7^)) + o(l), 

using that Tr(7„) A and Tr(7„) = Tr(7^)+Tr(72), which follows from xl+C = 1 
and the cyclicity of the trace. Because of lim„_j.oo £i"fm ««) = we deduce 

(7.19) lim £{-fi,al) = /(A^) and lim f (7,^,, a„) = /(A - A^), 

where we use that Tr(7,^J = J XnP-y„ ^ by Lemma |7.3[ and the the continuity 
ofAh^/(A). 

Therefore equation (|7.19p shows that both {(7;^,Q;^)} and {(7^,,a^)} are min- 
imizing sequences for /(A^) and /(A — A^), respectively. Since we already know 
that Tr(7^) — A^, we conclude that (jny'^n) ~^ (7^:Q!^) in from Corollary [ITT] 
By continuity of we also have lim„_j.oo ^^(7^, ck^) — £{j^,a^) = so that 

(7^,0;^) is indeed a minimizer for /(A^). This completes our proof of ()7.14p . 

Step 4: Binding Inequality and Conclusion. At this stage, we have proven 
that the energy of our minimizing sequence behaves like 

/(A) -/(Ai) + /(A-Ai) 

and we have seen that there exists a minimizer (7^, a^) for /(A^). Note that (7^, a^) 
is the weak limit of (7„, «„) and the strong limit of {jn,a^) = Xn{ln,cx.n)Xn in 
On the other hand, the sequence {(7^, a^)} is a minimizing sequence for /(A — A^). 
Hence one of the two following situations must occur: 

• either {(7^,0;^)} is relatively compact up to translations. In this case it 
converges (up to a subsequence) to a minimizer (7^,0^) for /(A — A^); 
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• or {(7,^,0^)} is not relatively compact up to translations. We may then 
apply the whole process again and we deduce that we have 

/(A - Ai) = /(A2) + /(A - Ai - A2), 

for some A^ > and that there exists a minimizer (7^, a^) G /C^2 for /(A^). 
We can summarize both cases by saying that 

(7.20) /(A) = /(Ai) + /(A^) + /(A - Ai - A^), 

for some A^,A^ > such that A^ + A^ ^ A and that I{X^) possess a minimizer 
(7^^, a'^) for k — 1,2. Next, we claim the following result. 

Lemma 7.4. (Binding Inequality). Assume that A^ and A^ are as above. Then 
one has 

(7.21) /(A^ + A^) </(A1)+/(A2). 
Using (|7.2ip and Lemma I4TT1 we deduce the strict inequality 

(7.22) /(A) ^ /(A^ + A2) + /(A - A^ - A^) < /(A^) + /(A^) + /(A - A^ - A^), 
which contradicts equation (|7.20l) . Hence the proof of Theorem [1] will be finished 



after we have proven Lemma 17.41 



(7.23) SiXRYXB., XRa'xB) ^ /(A^) + iov i = 1, 2, 



Proof of Lemma \ 7.4\ Consider two minimizers {j^,a^) and (7^,0^) for /(A^) and 
/(A^), respectively. We claim that 

for all i? > sufficiently large and some constant C > 0, where we recall the 
definition of the cutoff functions XR. s-nd Cr. from (jl.Sp . 

It suffices to show (|7.23p for i = 1, since the proof for i = 2 is analogous. As 
for the kinetic energy, the localization formula in Lemma [A. II yields Tr(T7^) ^ 
Tr(rxfl:7^Xi?) ~ C / . Next, we estimate the direct term as follows (using the 
symmetry in x and y): 

^ ff - xl{x)xliy))p^<x)p^iiy) ^^^^ 

ff p-yi{x)p^i{y) f C 
j j dxdy^CI p^i (x) dx , 

JJ{|x|^fl}xR3 \X — y\ J\x\^R. ^ 

for i? > sufficiently large. Here we used the Kato-Hardy inequality to conclude 
that * p^iH^oo ^ C\\^ p^i\\jji/2 as well as the decay estimate from Lemma 

15.31 Finally, we estimate the pairing term as follows: 

(l-XaW,Y«(i/)')l°'(i.ll)l' j^j^ 

JJ{\x\^R.}xR3 \x-y\ JJk3xr3 \x-y\ 

^ 7rTr(if Cfl/27'Cfl/2) ^cl f p^i (x) dx + Tr(rCfl/27'Cfl/2) ) • 

\J\x\^B/2 J 
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Here the first inequality follows from the symmetry |q;^(x,2/)| = |a^(2/,x)| and 
the support properties of XR- In the last line, we used Kato's inequality and 
{a^)*a^ ^ 7^. Putting all estimates together and using Lemma 15.31 again, we 
conclude that (|7.23l) holds. 

Next, we prove (|7.2ip . To this end, we consider the trial state 

(7.24) (7fl,ai^) := U{xRj'^XR,XRa^XR)U* + T5r^V{xrJ^XR,Xr(^'^XrW*T5r^, 

for some unit vector v and some rotations U^V €z 5*0(3). As XR &n.d tqr^xr have 
disjoint supports, it is easily seen that the above operator belongs to /C. Moreover, 
the number of particles of this trial state satisfies 

(7.25) Tr(7fl) = Tr(xfl,7'xfl) + Tr(xfl7'xfl) ^ + A^. 

Note the invariance by rotation and translation £(TyU{'j,a)U*T*) — £(7, a) for all 
?; S R'^ and all U £ 5*0(3). Since A i-> /(A) is decreasing by Lemma \4A\ we have 

/(A^+A^) ^ f f £ijR,aR)dUdV 

JS0{3} JSO{3) 

= £{xRi^XR,XRa^XR) + £{xri^xr,xr(^'^xr) 

I{x)g{y - 5Rv) 



y\ 



dx dy 



^/.o\ O k\ a 
</(Ai) + /(A^) + — - 



i?2 R ' 

which follows from Newton's theorem and (|7.23p . where we have defined 



fix) = / XRixYp-^i {U ■ x) dU, g{y) = / XRivY P-t^ {V ■ y) dV. 

JS0{3) JS0{3) 

This last estimate yields the desired result by taking R > sufficiently large. □ 
The proof of Theorem [T] is now complete. ■ 

8. Proof of Theorem H] 

The proof of Theorem |4] is a combination of ideas from the proof of Theorem 
[T] and of the usual geometrical methods for A^-body systems as was developed in 
m mm m ini [Sni- to our knowledge, they were applied for the first time 
to nonlinear models approximating the iV-body case (including the Hartree-Fock 
approximation) by Priesecke in |T3] . We shall only sketch the proof of Theorem S) 
In our opinion, it is much easier than the proof of Theorem [1] 

Step 1: First properties of HF energy. The first step consists in showing that 

(8.1) /"^(l) = 0, /"^(iV) < for all ^ 2 
and that 

(8.2) I^^{N) ==; /HF(7V -K) + I^^iK) for all K = L.A^ - I. 

The value of /^^(l) is obvious as for density matrices of rank one, the exchange 



and direct terms cancel. The rest is proved similarly as in Lemma 14.1 

The main part of the proof consists in showing that if 
(8.3) I^^{N) < mm{I^^{N - K) + I^^{K), K = L.A^ - 1). 
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holds true, then all minimizing sequences are compact and thus one gets the ex- 
istence of a minimizer for I^^{N). In the following we fix some N ^ 2, assume 
that (|8.3p holds, and consider a minimizing sequence {'jn} for I^^{N). Similarly 
to Lemma [42l we obtain that {(7^, 0)} is bounded in X. 

Step 2: No Vanishing. We claim that our sequence {7^} cannot vanish, similarly 
to Lemma 17.11 Indeed if {7^} would vanish, applying Lemma 17.21 we get that 
D{py^, p^^) — )• 0, hence we would obtain I^^{N) ^ which contradicts (|8.1I) . 
Therefore we may assume (translating 7„ if necessary), that 7„ ^ 7 with :— 
Tr(7) > 0. Also we have similarly to Lemma 17.31 that we can find a cut-off function 
on a ball of size i?„ such that V^x(^ri ^■)7nX(^n ^0 Vk^Vk strongly. 

Also, similarly to Corollarv l4.11 one shows that {7„} is compact when Tr(7) ~ N, 
hence it suffices to prove that = N. 

Step 3: A^-body Geometrical Methods. The next step is an easy adaptation of 
[15) , i.e. it consists in using a cluster decomposition like in the usual HVZ Theorem, 
based on the fact that 7„ arises from an A^-body wavefunction which is a Slater 
determinant. Indeed following jl5) . one arrives at an estimate of the form (up to 
extraction of a subsequence): 

N-l . . 

(8.4) y <^ (/"^(A) - /"^(A -K)- I^^iK) } lim ^ 0. 

K=l 

The nonnegative numbers are defined by 
/A^\ f f ^ ^ 

^K=[j^] / • • • / IlCn(2:j)^ Y[ Xn{xi)^\^„ixi,...,XN)\'^dxi---dxN, 
^ ^ •' •' i=l i=K+l 

where is the Slater determinant associated with 7„ and Xn — xi^n^') i^" 
troduced before, Cn = \/l — Xn> see [15]. Note that may be interpreted as the 
norm of the wavefunction associated with all the clusters for which K particles are 
sent to infinity while N — K particles stay close to zero. Note one has 

N N 

(8.5) y 1 and lim V (A - A:)zJ = A^ 

K=0 K=0 

The strict inequality (|8.3p and (|8.4p imply that we must have 

lim = for a: = I..N ~ 1. 

n— ^00 

This implies that we have 

lim / ••• / (n{xi)'^Xn{x2)'^\'i'n{xi,...,XN)fdxi---dxN ^0, 



n— )-C30 



lim 

n— J-oo 



^ — -v/z^'J^ — ■v/z^'I'^ 

V *n V^N^n 



0, 



where and are Hartree-Fock states defined as 
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Using that the operators 7*^ , 7^1 and 7^2 are projectors, we obtain 

Taking the trace of the previous equahty, we infer 

Um {z- - - hm (z]^ - (z]^)^) = 0. 

n— >-oo n— >-oo 

By (|8.5p . we have lini„_j.oo A^^o > 0, therefore we arrive at the following 

conclusion: 

hm = 1 and lim z^ = for = l..A^ 

n— >-oo n— >oo 

hence = iV. This proves that {(7n,0)} is indeed compact in X. 

Step 4: Proof of the binding inequality (|8.3p . We have seen that minimizing 
sequences are all compact when (|8.3p hold. It remains to prove (|8.3p . This is done 
like in il5, by induction on N ^ 2. As = and I^^i2) < = 2I^^{1), 

we deduce that /^^(2) possesses a minimizer. Next the strict inequality /^^(3) < 
/H^(2) + = /"''(2) is proved as in [15, Thm. 3], using that I^'^{2) possesses 

a minimizer. The rest follows by induction, as soon of course as the energy stays 
bounded from below and coercive, i. e. when N < N^'^{k) holds. 

This completes our sketch of the proof of Theorem S) ■ 



Appendix A. Localization of Kinetic Energy 

In the section we collect some localization estimates for the pseudo-relativistic 
kinetic energy operator \/— A + m^. In fact, similar localization estimates - often 
referred to as IMS-type localization formulas ~ can be found, for instance, in the 
stability analysis of relativistic matter [33) . However, the following IMS-type esti- 
mates for \/~-A + m? turn out to be suitable when studying lack of compactness 
of minimizing sequences. 

Lemma A.l. Let {xi}iG/ be a smooth partition of unity with = 1 

such that each € L^(M'^) with p G {2d, 00]. Then 

where K = v^A + tt? with m > 0. 

The proof of Lemma [A. II readily follows from the next result. 

Lemma A.2. Suppose x e C°°(M'') such that Vx e LP{W^) with p e {2d, 00], and 
let K — A -f- with m > 0. Then the following formula holds: 

(A.2) -{Kx^^X^K)=xKx--l J^.\^X?J^,V-sds^T 
where is a nonnegative bounded operator. Moreover, we have 

r°° ^ 1 _ 

^C\\Vx\\U 



(A.3) / j^\^x\^j^,Vrsds 



6p/2(i^(K'')) 
2 



for some constant C independent of x (^^t depending on m, p and d). 
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Proof. By functional calculus, we have the representation formula 



K = 



1 



ds 



s + K^ Vs' 



(A.5) 
which leads to 

1 Z""" 1 1 
(A.6) [^'^] = -| J^i^'^B]-^Vrsds. 

Using (jA.6p . a calculation shows 



i {Kx^ + X^K) = xKx +l[[K,x], X] 

1 



1 



= xKx 



1 



1 



27r 7o s + K"^ 



[[K\xlx, 



s + K^ 
1 



X 



with 
(A.7) 



s + K^ 
1 



-y/s ds 
y/s ds + L^, 



if 2 I 



if2 



■\/s ds. 



Since [[if^^ -j^]^ -^J = — 2|VxP, this yields (jA.2|) . Furthermore, we note that 



(A.8) 



= - 



A{syA{s)^/sds, A{s) 



{s + K^y/^ 

which shows that is nonnegative. 

We now show that all the above commutators are bounded (which will also 
validate our somewhat formal calculation). Recall the Kato-Seiler-Simon inequality 
(see [31] and [Ml Thm 4.1]) 

(A.9) ||/(-zV,).g(a:)||e^(^.(R.)) ^ {2n)-''/P \\f\\^,^^.^ , for p S [2, oo]. 



This yields, choosing some d/{2p) < q < 1, the following bound. For all s ^ 0, 
1 „ 1 



(A.IO) 



s + K^ 



s + K^ 



ep/2(i^(R'')) 
C 



_ ||('s + f2^-9||^ 



2 



Furthermore, a scaling argument shows that 



2^-91 



Hence we obtain 



(A.ll) 



1 



S + if2 



X2 



gd/(2p)-g 11(1 + ^2-) 



2^-9|| 



i ds 



^C\\\7x\\l.(m 



oo l/2+d/p-2g 



■ ds. 



where the last integral is convergent under our assumptions on d and p. The 
argument is exactly the same for L^, noting that [K'^,x\ — ~i{P ■ Vx + Vx ■ P) 
with P = -iV and using that \\P{s + + \P\'^y^''^\ ^1. □ 
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Appendix B. Proof of Proposition 12.1 
We consider the following minimization problem 

(B.l) /o,.(A)= inf \TriV^j)-!^D{p„p^) 

^2^^^^'^ 2 77 \x-y\ 

By rescaling j{x,y) i— > j3'^^{l3x, j3y) and a{x,y) i— >■ P^a{Px, Py) with /? > 0, we 
deduce that either /o,k(A) — or — oo. Following the proof of Lemma BTTl one also 
sees that /q is non-increasing. Hence there is a uniquely defined A^^^(k) such 
that /o,k(A) = if and only if A ^ A^^^(k). Next, we notice that 



A — TO ^ V^ — A + TO^ — TO ^ %/ — A 

hence /o,k(A) — toA < /(A) ^ /o,k(A). This shows that /(A) is bounded below if and 
only if A ^ A^^^(k). Furthermore, it is clear that A^^^(-) is nonincreasing, since 
for all (7, a) G JC, 

-D{p,,p,) + Ex(7) - jj ^^dxdy ^ 

by (|2.1ip . The continuity of A^^^(k) is easy to verify. 

Next, we prove that A^^^(k) > and we derive the asymptotic behavior of 
A^^^(k) as k — ^ stated in Proposition 12.11 To this end, we introduce another 
homogeneous minimization problem 

(B.2) /5l'(A) = , inf^ {Tr(y=A7)-fD(p^,p,)}. 

which, of course, yields a corresponding We claim that A^^^(k) ^k-s-o 

y^^^ifi). Indeed, by the nonnegativity of the exchange term and (|2.16p . 

(B.3) /o..(A)>(l-^)/S:^(A). 

Hence 

(B.4) A"™(k) ^ A'''^'^ 



1 - K7r/4^ 

Now we recall that A''°'^(k') > for every k' > 0. Indeed we have 

(B.5) / / ^iM^dxd, ^ C\\p,fl:t\\p,\\% ^ CA^/^Tr(^/3A^) 
JJr^xr^ \x-y\ ' 

where we have used the Hardy-Littlewood-Sobolev inequality and the well-known 
[TT] semi-classical estimate Tr(-y/— A7) ^ ||p-^||^{^3 for fermionic density matrices 

< 7 ^ 1. This shows that A'<=^(k') ^ {2/{Cn')f^^ > 0. Hence we deduce from 
(|Rl| and < K < 4/7r that A"p^(k) > 0. 

Next we derive an upper bound for A^^^(k). We have for any (7,0) G ICx, 

(B.6) /o,«(A) ^ f (7, 0) = Tr(y^7) - ^D{p^,p^) + ^Ex(7). 

Now we have similarly to (|2.16p 

(B.7) // ^^^dxdy ^ |Tr(y3A,.) ^ ^TriV^,). 
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Hence we obtain 

(B.8) /o,4A)5;(l + ^)/5r^(A), 

and we thus get the upper bound 

(B.9) A«™(«)<A-<if— ^ 

The last step consists in showing that V^'^{k) behaves as stated when k — >■ 0, 
which will imply the result for A^^^(k) by (|B.4p and (|B.9p . This was indeed essen- 
tially done by Lieb and Yau, see Corollary 1 in [32] ■ • 



Appendix C. Proof of Proposition 14.2 



We want to prove that G(A) = fiX, where /3 is the first eigenvalue of the operator 
K—k/2\x\ acting on iQjj^(R'^) when q = 1 and L^(]R^) when g > 2. For the reader's 
convenience, we recall that K — \/~A + m"^ . 

We start by deriving a lower bound. Using ()2.15p . we infer Ti{K^) ^ Tr(ifaQ:*) 
and thus 

K ^ \a{x,y)\'^ 

- — a, 2 
Tr(Q*a)^A 

The right-hand side of the previous inequality can also be written, using the anti- 
symmetry of q;, as follows. 



(C.l) 



G(A) ^ inf <^ Tv{Kaa 



-dx dy 



(C.2) 



inf 

3^(a:,y)^ = -a(2;,y), 



2 X- 



II "II 



This quantity is just A times the lowest eigenvalue of the two-body operator 
^^•^^ ~~2 2^ 

when it is restricted to the fermionic space L^(M'^,C'^) A L^(M'^,C'^). If 9 = 1, 
its lowest eigenvalue equals the first eigenvalue of the same operator acting on 
L^(]Ri^,M) A L2(M^^]R) (triplet state), whereas when g ^ 2, it coincides with the 
first eigenvalue of the same operator acting on L^(R^,R) (gjg L^(M^,R) (singlet 
state). Removing the center of mass as explained in [27], one sees that this lowest 
eigenvalue is given by 



inf inf (Tj--, 




2\u\ 



where u ~ x — y and i is the Fourier variable associated with the space variable 
V ^ {x + y)/2. We have to take S) = Ll^^{m!^) when q = 1 and 55 = ^^(rS) ^^^^^ 



q ^ 2. The infimum above is attained for i — 0, and hence we obtain the desired 
lower bound: 



(C.4) 



G(A) ^ XMay, v/-A„ + - — := /3X. 



2\u\ 



Note that we can in fact deduce that the infimum on the right side is attained. 
This follows from the min-max principle and the operator inequality V— A + — 
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2\x\ 



pr, where the right side is the Schrodinger operator for the 



2\x\ 



nonrelativistic hyrdogen atom. 

Next, we prove the upper bound in a similar way as done for the Miiher functional 
in [M]. Since K—k/2\x\ is a real operator, there is a real- valued eigenfunction f £ S) 
such that 



(C.5) 



K 

2bi 



/-/?/• 



Using the method in TU], one can show that / € C°°(]R'^ \ {0}) a-nd that it decays 
exponentially at infinity when ^ k < A/tt. Note that / e D{\/—A + m? — 



K/(2|a:;|)) C i7^/^(R^) but this domain is only known to be _ff^(M^) when k ^ 1. 
(We will address this regularity issue below again). 

Next, we fix a nonnegative function x G C^(R"^) such that J x'' = A, and we 
denote xl{x) = L^'^/'^xi^l L) for L > given. We define 



(C.6) 

where S 



1 if g: 



0iL{x, y) = XL{x)f[x - y)xL{y)^, 
1 and, when q ^ 2, 



1 

V2 



/O 1 

-1 



V 



0/ 



As aL e L'^i 



^ C X C), we may denote by aL the Hilbert-Schmidt operator 



whose kernel is aL{x, y). Then aj^ = —aL by the choice of /. Note that 



(C.7) 



\\XL\\l^=OiL-'^') 



and hence — >■ in the operator norm. Next, we define jl as the unique nonneg- 
ative trace-class operator solving the following equation for L > large enough: 



(C.8) 



7-l(1 - 1l) = aLa*i 



Note that 7^ < 2aLa'L if L > is sufficiently large. In particular, we have ||7l|| — 
0{L~^). In analogy to 14 , it can be checked using ^ x^ = ^ that Tr(aia^) — 
A -t- 0{L-'^). Therefore, we have Tr(7L) = A + 0{L~'^), by (jOSl) . and the fact that 
Tr(7i) = 0(L-3). 

Next, we show that {'Jl, c(l) G A^- To see this, we recall from W and Prop. 
2] that there exists an orthonormal basis {(pi) of L^(R^,C^) such that 



(C.9) 



aL 



On the other hand, the operator a^a^ is just 1/2 times the one-body density matrix 
of the two-body wavefunction a^, so that 



(C.IO) 



2i-l I 



\^2i){(P2t\) , 



i>l 



(C.ll) 



iL^^li {\f2z-l){'fi2i-l\ + \'fi2z){(P2i\) , 
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where for all i ^ 1, — 7i) — (ci)^/2. Next, we have to verify that 



(C.12) 



7L 



1 - 7l 



s; 1. 



As aL^2i-i = —-^V2i and aL^2i = :^</'2i-i with respect to the basis (v^i); this 
amounts to studying the following matrix 



(C.13) 



M = 



( l^ c,/V2^ 

7, -c,/V2 

-c,;/\/2 1-7, 
\c^|^^2 1-7,/ 

which is easily shown to satisfy ^ A/ ^ 1 (it is indeed a projection). Thus we 
have proven that (7l, a^) S /C for L > sufficiently large. Next, we estimate 



(C.14) 



K f f |aL(x,y)p 



\x - y\ 



■dx dy + Tr{K-fl). 



We start by discussing the kinetic energy. Using (jA.2p . we obtain 
Tr{KaLal) < (2^)3/2Tr [xlK{p)f{p)xlf{p)] 

(2^)3/2 



; (isTr 



xi 



Since the eigenfunction / satisfies ()C.9p . we obtain 

\aL{x,y)\'^ 
\x - y\ 



{2nf/'TT[xlK{p)f{p)xlf{p))=^ 



Notice the above expressions are well-defined even when / only belongs to H^^^{ 
When 4/7r > K > 1, they can be justified by approximating the eigenfunction / in 
7^1/2 ^]^3 -J gjjiooth functions and passing to the limit. On the other hand, we 
have 



fip) 



K(pf 



fip) 



K{pf 



-|Vxi 



K{pf 



fip) 



K{pf 



-xl 



II 



WxlW 



L4 



62 



m + s 771 + 5 ~ (771, + s)2 

Moreover, for L > sufficiently large, 

(C.15) Tr(i^7i) s; 2Tv{K{aLalf) ^ 2 Tr(KaLal) = 0(£-'). 

In summary, this shows Gi^LjCtL) ^ f3X + 0{L~^). Passing to the limit L — > 00, 
we get the upper bound G{X) ^ /3A. Combined with the lower bound, this finally 
yields the desired equality G(A) — (3X. 

The estimate for /(A) is obtained using the same state (7^,0^). We have 

^f(^rf.rf,^fTr(if7i) = 0(X-). 
\x-y\ 2 

which we proved above in (jC.lSp . Moreover note that p^^ ^ Poia' by (jCSp . and 
this is seen to imply 
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as can be deduced from pL4t p. 052517]. Therefore, by choosing L > sufficiently 
large, we obtain 

/(A) ^ G(A) -mX- ^D{x\ x') + o{L-') < G(A) - ™A, 

as desired. ■ 
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